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ABSTRACT 

Although photometric and spectroscopic surveys with the Spitzer Space Telescope increased re- 
markably the number of well studied debris disks around A-type and Sun-like stars, detailed analyzes 
of debris disks around F-type stars remained less frequent. Using the MIPS camera and the IRS 
spectrograph we searched for debris dust around 82 F-type stars with Spitzer. We found 27 stars that 
harbor debris disks, nine of which are new discoveries. The dust distribution around two of our stars, 
HD 50571 and HD 170773, was found to be marginally extended on the 70/im MIPS images. Combin- 
ing the MIPS and IRS measurements with additional infrared and submillimeter data, we achieved 
excellent spectral coverage for most of our debris systems. We have modeled the excess emission of 22 
debris disks using a single temperature dust ring model and of 5 debris systems with two-temperature 
models. The latter systems may contain two dust rings around the star. In accordance with the 
expected trends, the fractional luminosity of the disks declines with time, exhibiting a decay rate 
consistent with the range of model predictions. We found the distribution of radial dust distances as a 
function of age to be consistent with the predictions of both the self stirred and the planetary stirred 
disk evolution models. A more comprehensive investigation of the evolution of debris disks around 
F-type stars, partly based on the presented data set, will be the subject of an upcoming paper. 
Subject headings: circumstellar matter — infrared: stars 



I. INTRODUCTION 

Nearly all young stars harbor circumstellar disks, 
which initially serve as a reservoir for mass accretion, 
and later can become the birthplace of planetary systems. 
During this latter process, the originally submicron-sized 
dust grains start growing, and their aggregation is be- 
liey cd to lead to km-size planetesimals (for a review, 
see lApai fc Laurett a 2010, and references therein). Non- 
destructive collisions between planetesimals result in the 
formation of subsequently larger bodies. These events 
happen first in the inner disk due to the shorter coUi- 
sional timescales, then the process propagates outwards 
(jKenvon &: Bromleyl l2004ah . The newly formed Pluto- 
sized protoplanets stir up the motion of leftover smaller 
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bodies in their vicinity, initializing a collisional cascade. 
As they become more energetic, collisions result in the 
erosion of planetesimals and the production of small 
dust grains. An optically thin debris disk is formed, in 
which the second generation, short-lived dust grains are 
continuously replenis hed by collisions and/or evapora- 
tion of planetesimals (jBackman fc Parescelll993t iWvatti 
120081 ). This self-stirring mechanism is not the sole fea- 
sible way to incite destructive collisions between minor 
bodies. Giant planets, formed previously in the pri- 
mordial disk, or stellar companions can also dynami- 
cally excite the motion of planetesimals via their secu- 
lar perturbation, even at a significant distance from the 
planetesimal disk. Thus these large bodies can also ini- 
tiate and gov ern the for mation and evolution of a de- 
bris disk (Mu still fc Wva tt 2009), providing an alterna- 
tive stirring mechanism. In a debris disk, mutual col- 
lisions grind down planetesimals to small dust grains 
that are then ejected by radiation pressure, or in more 
tenuous disks re moved b y the P oynting-Robertson drag 
llDominik fc De cin 20031: IWvatti [2005). This process is 
accompanied by the depletion of the reservoir planetesi- 
mal belt an d eventually leads to the decline of the debris 
production (jWvatt et al.ll2007al: iLdhne et al.ll2008D . 

Due to the strong link between the debris dust and 
the unseen planetesimals, the investigations of the small- 
est particles of debris systems can lead to a better un- 
derstanding of the formation and evolution of planetes- 
imal belts and, eventually, the formation and evolution 
of planetary systems. The observational verification of 
the different aspects of planetesimal formation and evo- 
lutionary model predictions requires a detailed study of 
the incidence of stars with infrared (IR) emission due to 
debris dust and investigating the change of debris disk 
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properties (e.g. radius of the dust ring, fractional lu- 
minosity) with age. The ideal way would be to resolve 
and observe many debris disks in scattered light or in 
thermal emission from optical to millimeter wavelengths 
with good wavelength coverage. In reality, however, the 
number of resolved disks is very limited and the spec- 
tral energy distribution (SED) of the dust emission was 
measured for most debris systems only in a few infrared 
bands. The fundamental parameters of the disks have 
to be estimated from these sparsely sampled SEDs. The 
interpretation of SEDs is ambiguous (e.g. considering 
the radial location of the dust) but by handling a debris 
disk sample as an ensemble, one can obtain a meaning- 
ful picture about the basic characteristics of the parent 
planetesimal belt(s) and about the evolutionary trends. 

The current theoreti cal models dealing with th e build 
up of planetesimals (Kenvon fc BromlevI l2004aL |2008| ) 
and with the stea dy-state collisional evo l ution of the 
planetesimal belts (Dominik fc Declnll2003l : iWvatt et al.l 
I2007at iLohne et al.ll2008D predict how the fundamental 
properties of debris disks evolve with time. At a spe- 
cific radius, the peak of the dust emission is believed 
to coincide with the formation of 1000-2000 km sized 
planetesimals. After this stage - parallel with the de- 
pletion of planetesimals - the dust emission decreases 
with time. The evolution of the disk can be traced both 
in the variation of incidence of disks with time and in 
the evolution of the brightness of dust emission. The 
different models predict that the dust fractional lumi- 
nosity, the ratio of the energy radiated by the dust to 
the stellar luminosity, varies with time as t~'^, where 
n = 0.3 — 1 in disks where collisions are t he dominant 
remov al process fo ominik & Dccin 200 3t IWvatt et ahl 
l2007atlL5h^e et aLir2008l : lKenvon fc Bromlevll2008[ ). The 
unique sensitivity of the Spitzer Space Telescope in the 
MIPS 24/im band allowed the detection of stellar photo- 
spheres and a sma ll amount of exces s for a large number 
of field stars (e.g. iRieke et al.|[2005t lMever et al.1 [20081) 
and even for relatively distant open cluster m embers 
(e.g.l ^ung ct al. 2004; Gor lova et al . 2006; Sicgle r et aP 
[2007noilrrie et al.l 120081 : iBalog et al.l 120091 1. The lat- 
ter observations enabled the study of the evolution of 
warm dust around well-dated sampl e stars. Investigat - 
ing early (late B- and A-) type stars IRieke et al.l ()2005D 
demonstrated a decline of debris disks with age: older 
stars show excess emission less frequently and wi th lower 
fractional excess than the younger ones. Siegl er et al.l 
■|2()07 ; found similar evolutionary trends for debris disks 



encircling later type stars (F,G,K). Based on obser- 
vations of more than 300 Su n -like s tars with spectral 
type of F5-K3, iMever et al] ()2008D argued that the 
24/zm excess fraction for this sample is roughly constant 
for ages <300Myr and declines thereafter ( see also in 
iCarpenter et al.H2009aD. Re centlv. iGas par et~all ()2009l ) 
and lCarpenter et al.l ( 2009b[ ) gave a summary of the evo- 
lution of 24/^m excesses around B7-K5 type stars. Con- 
firming the previous results they concluded that both 
the incidence of 24/im excess and the excess luminos- 
ity monotonically decrease with time at ages >20Myr. 
Utilizing the observations of different infrared space mis- 
sions {IRAS, ISO, Spitzer), the predicted evolutionary 
t rend in the fractional luminosities was also established 
Decin et al 
Rhec et al.. , 



2003"; 'Su et a ll l200l IWvatt et al.l I2007bl: 
2007; .Carpenter et al.l l2009aD . 
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tended planetesimal disk both the stirring by Pluto-sized 
planetesimals that were born in the same belt, and the 
dynamical excitation by secular perturbation of distinct 
giant planets is thought to be accompanied by the out- 
ward p ropagation of the dust production site with time 
(iKenvon fc Bromley 200.8.; Mus tiU fc Wvat t 2009) . The 
observational evidence for such a delayed initiation of 
the collisional cascade as the function of radial loca- 
tion is not yet conclusive. Some surveys did not report 
any trend in the e volution of the radius with age (e.g. 
Naiita fc Williaini 2005), wh ile studying debri s disks 
around B- and A-type stars iRhee et aJj (|2007l ) found 
some evidence that the radius of dust belts is increas- 
ing with stellar age. 

Thanks to the recent photometric and spectroscopic 
surveys with the Spitzer Space Telescope, the num- 
ber of debris disks with detailed spectral energy dis- 
tribution at mid- an d far-IR wa velengths has been in- 
creased significantlv (iChen et al.li 2006; Riokc ct al. 200{| 
Su et al.1'2006'; 'Car penter eiTaiHoOS ; .RebuU et al.„2008l 
Trilling ct al. 2008). This improvement is especially re- 
markable for disks around A-type and Sun-like stars (late 
F, G, and K-type stars). The comparison of these data 
with the predictions of quasi-steady state evolutionary 
models showed that most observed trends for A-type 
and Sun- like stars can be rep r oduced adeauatelv (iWvatti 
l2008t ICarpenter et al.ll2009al) . IKennedv fc WvattI (|2010[ ) 
confronted the Spitzer observations of A-type stars with 
an analytic model that also take into account the effects 
of the self-stirring on the disk evolution. Utilizing this 
model they were able to reproduce the observed trends 
and they obtained rough estimates for some initial pa- 
rameters (e.g. average mass) of disks around A-type 
stars. It was also concluded that debris disks are nar- 
row belts rather than extended disks. According to the 
models, F-type stars are expected to be an intermedi- 
ate type between the A-type and Sun-like stars in terms 
of debris disk evolution as well: 1) their disks are pre- 
dicted to evolve faster than those around main-sequence 
stars of later types (in disks with identical surface den- 
sity distribution, the timescale of planetesimal formation 

processes are thought to be proportional to M* ' ); 2) F- 
type stars live much longer than A-type stars (the main- 
sequence lifetime of an 1.4 M0 F5-type star is 3 times 
longer than the main-sequence lifetime of a 2.0 M0 A5- 
type star), making it possible to follow disk evolution 
for a significantly longer time. Up to now the number 
of detailed studies of debris disks around F-type stars is 
modest compared to the A-type and Sun-like samples, 
preventing us from understanding the evolutionary as- 
pects. 

In this paper we present the results of a large survey 
with the Spitzer Space Telescope that focuses on debris 
disks around F-type stars. Our main goals are to 1) 
significantly increase the number of debris disks with de- 
tailed SED around F-type stars; 2) investigate the varia- 
tions of fundamental properties of the disks, and compare 
the observed trends with the predicted ones; 3) compare 
the evolutionary trends obtained for disks around A-, F-, 
and G/K-type stars. In the present paper we review the 
target selection (Sect. [2]), observations and data reduc- 
tion aspects of the F-stars program (Sect. |3]). We iden- 
tify stars with infrared excess, model their SED and esti- 
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mate the fundamental properties of the observed debris 
disks (Sect. [J]). Using the derived parameters we inves- 
tigate the diversity of the fundamental disk properties, 
and compare the observed trends with the predictions 
(Sect.©. 

Four new warm disks - discovered in the framework of 
this program - hav e already been analyzed and published 
()Mo6r et al.ll2Q09[ ). The evolutionary aspects of the cur- 
rent data set - supplemented by the recently discovered 
four warm debris systems, as well as additional debris 
disks around F-type stars observed by Spitzcr from the 
literature - will be further analyzed in an upcoming pa- 
per (Moor et al. 2010a, in prep.). 

2. SAMPLE SELECTION 

Our primary intention was to study the variations of 
the disk properties, thus in the sample selection we fo- 
cused on stars where previous IR observations hinted on 
the existence of excess emission. Although the formation 
of planetesimals may last as much as hundred million 
years in an extended disk, the most active period of this 
process is restricted to the first few tens of millions years 
pCenyon & Bromley 2008). The possible outward prop- 
agation of the planetesimal formation can be verified in 
this period. Since the age of nearby moving groups over- 
laps well with this period, they are favorable, nearby, and 
well-dated places for investigations of the debris disk evo- 
lution process. Thus, the above mentioned sample was 
supplemented by several F-type members of the nearby 
young kinematic groups. For these young stars we did 
not require a priori information about the presence of 
emission in excess to the stellar photosphere. Because of 
the selection method the sample is inherently biased with 
respect to the presence of disks, therefore it can not be 
used to study the incidence of debris disks around F-type 
stars. 

With the aim of constructing a list of F-type main- 
sequence (and some subgiant) stars, where earlier ob- 
servations indicated the presence of mid- and/or far- 
infrared excesses, we carried out a systematic search us- 
ing the data of the IRAS and Infrared Space Observatory 
satellites. For the selection of IRAS-based candidates, 
we collected all sources from the IRAS Faint S ource 
Survey Catalogue (IRAS FSC; IMoshir et all [19891) and 
the IRAS Serendipitiou s Survey Catalogue (IRAS SSC; 
iKleinmann et al.l [1986^ having at least moderate flux 
quality at 25 and 60/xm. The positions of the IRAS 
sources were cross correlated with the entries of the Hip- 
parcos Catalogue (ES A 1997) and the T ycho-2 Spec- 
tral Type Catalogue (jWright et al.l |2003[) and we se- 
lected only those objects where the positional match 
was within 30". Giant stars were omitted from the 
sample on the basis of the spectral information or the 
absolute magnitude. For the id entification of stars 
with e xcess we used the method of IMannings fc Barlowl 
(|1998D . First, we computed the ratios of measured 25^m 
and 60/xm flux densities to the measured flux at I2fim 
(^12/25 = Itf > -R12/6O = ItJ) as '^ell as the error of 

the ratios = \J i^^f + ^ '^^12/60 = 

TTo \/(^)' + (^)') ■ Then, applying a Kurucz model 
atmosphere of a typical F-type star (a model with ef- 
fective temperature of 6500 K, log (7=4.25, solar metal- 



licity) we derived the expected photospheric flux ratios 
for the specific IRAS bands (^12/25' -^12/60)- calcu- 
lated the significance of the differences between the ob- 
served and expected ratios as 512/25 = ^^s^Ri2/2^^^^ ^^'^ 

'5'i2/60 = ^^^mi2^eo^'^° ' ^i'^ahy, we selected those objects 
where 6*12/25 < —2.5 or/and <S'i2/6o < —3.0. Due to the 
low spatial resolution of IRAS, many of the positional 
coincide nces between a star and an infrared source could 
be false ()Mo6r et al.ll2006HRhee et al.ll2007[ ). In order to 
reject objects that are possibly affected by source con- 
fusion we applied the same criteria as those adopted in 
Sect 2.3 in lMoor et all £2006). 

Observations with ISO confirmed the existence of IR 
excess and the positional agreement between the opti- 
cal and IR source for several of our IRAS-based can- 
didates. Moreover, the compiled list w as further sup - 
pleniented by five ISO- based discoveries (jSpangler et al.l 
[2OOII: IDecin eFall [200I . 

The list of F-type stars belonging to different nearby 
young stellar kinematic groups (e.g. Tucana-Horologium 
associatio n, AB Dor moving group) wa s adopted from the 
catalog of iZuckerman fc Song! ()2004b '). This initial list 
was amended by adding several new unpublished moving 
group members (Moor et al., 2010b, in prep.). 

After merging the lists, we excluded all those objects 
that were reserved by other Spitzer programs. Among 
the finally selected 82 candidates, 27 were based on hints 
for excess emission at 25/im (Warm Disk Candidates, 
hereafter WDCs), 46 are suspected to exhibit IR excess 
at longer wavelengths (Cold Disk Candidates, hereafter 
CDCs) and 9 were selected because of their kinematic 
group membership (Moving Group Members, MGMs). 
Note that many stars selected on the basis of their sus- 
pected excess are also young kinematic group members. 
The basic properties of the sample stars as well as the 
reason for their selection are listed in Table [T] 

2.1. Basic properties of the selected objects 

In order to estimate some basic properties of our stars 
and to provide photospheric fiux predictions at relevant 
mid- and far-IR wavelengths, we modeled the stellar 
phot osphere by fitting a n ATLAS9 atmosphere model 
( Cast elli fc Kuruc3l2003[ ) to the optical and near-IR ob- 
servations (Hipparcos, Tycho2, 2MASS). The surface 
gravity value was fixed during the fitting procedure. If 
there was no indication that a star has already left the 
main-sequence - based on its position in the HR diagram 
and/or spectroscopic information - we adopted a value 
of log g =4.25 corresponding to the main-sequence stage. 
For evolved objects the logg values were either taken 
from literature or computed from available data (see Ta- 
ble [1] for details). The metallicity data were also col- 
lected from the literature. In the cases where more than 
one [Fe/H] estimates were available for a star, we used 
the average. If no metallicity information was found for 
a specific star, we adopted solar metallicity. Most of our 
stars are located i nside the Local Bubble , where the mean 
extinction is low ([Lallement et al.ll2003D . Thus, for stars 
within 80 pc the visual extinction was neglected (Ay set 
to 0.0). For more distant objects and for stars without 
reliable distance information, the Ay value was a free pa- 
rameter in fitting the photosphere. Our sample contains 
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23 multiple systems (see Table [T]) . For close binary sys- 
tems with separation <5", we used the combined pho- 
tometry of the components, except for HD 122510 and 
HD 199391, where good quality photometric data were 
available separately for each component. In the latter 
cases the fitting was performed for the different compo- 
nents separately and then the obtained photosperic mod- 
els were co-added. The fixed (logg, [Fe/H]) and fitted 
(TcSyAy) parameters are quoted in Table [TJ 

3. OBSERVATIONS AND DATA REDUCTION 

We performed observations usin g the Multiband I mag- 

l|200l and 
on 



ing Photometer for Spitzer fMIPS .IRieke et al 



the Infrared Spectrograph f IRS. iHouck et al 



200l 



the Spitzer Space Telescope f Werner et al.l 120 04^ in the 
framework of the programs FID: 3401, 20707, and 40566 
(FI: F. Abraham). The MIFS photometry of HD 35114 
was taken from the Spitzer Archive (FID: 3600, FI: I. 
Song). The MIFS images were obtained in photomet- 
ric imaging mode (default scale, small- field size). All of 
our objects were measured at 24/im and 70/im, CDCs 
and most MGMs (in total 52 objects) were observed at 
160/im as well. In the 24/xm band we performed 2 obser- 
vation cycles with integration time of 3 s providing typ- 
ically 28 DCE frames per source. At TOfim and 160/xm 
2-10 cycles were made, the integration time was set to 
3 s or 10 s depending on the actual source. We obtained 
low-resolution spectra, using standard IRS staring mode 
observations with the Short-Low and Long-Low modules, 
which are sufficient to outline the continuum and the 
broad spectral features. Whenever it was possible, peak- 
up observations in high accuracy IRS or FCRS mode 
were used to place the star at the center of the slit with a 
radial pointing uncertainty of 0.4" or 0.14", respectively. 
For most WDCs/MGMs and for some CDCs the full 
wavelength coverage was requested (all four IRS mod- 
ules were used), while in the remaining cases the short- 
est wavelength module (SL2, 5.2-7.4 /.tm) was dropped. 
In the case of HD 122106 and HD 125451, where no high 
accuracy peak-up observations were available, only the 
longest wavelength modules (LL2, LLl) were utilized 
(at long wavelengths the slit is larger and therefore the 
probability of loosing flux due to inaccurate position is 
expected to be lower). For the two brightest sources 
(HD 15745, HD 170773) low resolution MIFS SED spec- 
tra (55-95/im; A/AA ~ 15-25 ) were also obtained. Six 
observing cycles were performed in both cases. Each cy- 
cle provided six pairs of 10 s long on- and off-source ex- 
posures. The on and off positions were separated by 2'. 

3.1. MIPS 
3.1.1. Image processing 

The processing of MIFS data started from the Spitzer 
Science Center (SSC) basic calibrated data (BCD) pro- 
duced by the pipeline version S16.1. Additional pro- 
cessing steps, including a fiat field correction and a 
background matching at 24 /xm, were performed us- 
ing MOFEX (Mosaicking and Foint source Extraction, 
IMakovoz fc Marleau|[200l . At 70 fim column mean sub- 
traction and time filtering were app lied for the images fo l- 
lowing the steps recommended bv iGordon et al.l ()2007[ ). 
Finally, the improved BCD data were coadded and cor- 
rected for array distortions. Bad data fiagged in the BCD 



mask files, as well as permanently damaged pixels flagged 
in the static pixel mask file, were discarded during the 
data combination. Outlier pixels were rejected using a 5, 
3 and 3 a clipping thresholds at 24/im, 70/im and 160/im, 
respectively. Output mosaics had pixels with sizes of 2'.'5, 
4" and 8" at 24/im, 70/im, and 160/im, respectively. The 
data reduction of the MIFS SED observations was also 
started with the BCD images (pipeline version S16.1) 
and the MOFEX was utilized to perform the necessary 
processing steps (combination of data, background re- 
moval, application of dispersion solution) and to compile 
the final image with pixel size of 9.8". 

3.1.2. Photometry 

We used a modified version of the IDLFHOT routines 
to detect sources and to obtain photometry on the mosaic 
images. 

MIPS 24. — All of our sources were detected with good 
signal-to- noise (S/N>20) on the final 24:iim images. Fig- 
ure [1] shows the positional offsets between the obtained 
centroids of the identified point sources and the 2MASS 
position (corrected for the proper motion between the 
epochs of the two observations). None of the sources 
show more than 3cr deviations from the averaged posi- 
tion. Apart from HD 38905 - offset of 1.42" from its 
2MASS position - the angular offsets between the dif- 
ferent positions are within the la uncertainty (~1.4") 
of the pointing reconstruction at 24/im (see MIFS In- 
strument Handbook). Thus, we verify that the observed 
24/im emission is emerging from our targets. 

Aperture photometry was performed to estimate the 
flux densities of the targets. The aperture was placed 
at the measured centroid of the source. The aperture 
radius was set to 13" and the background was com- 
puted in a sky annulus between 20" and 32". In the 
course of the sky estimates, we used an iterative sigma- 
clipping method, where the clipping threshold was set 
to 3(7. An aperture correction factor appropria te for 
a stell ar photosphere was taken from Engclbracht et all 
(l2007l) . In some special cases we had to deviate from 
this standard method. For four sources (HD 36248, 
HD 143840, HD 185053, HD 218980) that are surrounded 
by nebulosity at 24/im or/and at 70/im, a smaller aper- 
ture with radius of 3'.'5 was used. In cases where a 
bright nearby source contaminated the aperture photom- 
etry (HD 34739, HD 38905, HD 145371, HD 184169), we 
performed FSF ph otometry. The model FS F was con- 
structed following lEngelbracht et al.l ()2007D . The un- 
certainties of the photometry were computed by adding 
quadratically the internal error and the absolute calibra- 
tion uncertainty of 4% (MIFS Data Handbook). Due to 
the fact that our targets were measured with high signal- 
to-noise in this band, the uncertainty of the absolute cal- 
ibration dominates the error budget. The resulting 24/im 
photometry is presented in Tabled 

Figure [2] shows the distribution of flux ratios - mea- 
sured flux densities {F24) relative to the predicted fluxes 
{P24) - at 24/im as the function of the predicted photo- 
spheric fluxes. The histogram of the F24/P24 ratios shows 
a peak at around unity with a mean of 1.00 and disper- 
sion of 0.038. The obtained dispersion is consistent with 
the quoted absolute calibration uncertainty. Stars with 
F24/P24 >1.12 have excess emission at this wavelength. 
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MIPS 70. — At 70/im 30 among our 82 targets were de- 
tected at > 4(7 level. 

Aperture photometry with an aperture of 16" and sky 
annulus between 39" and 65" in radius was performed 
for all 82 targets. The sky level was estimated using the 
3a clipped mean. For the detected sources the aperture 
was placed around the derived centroid of the object, 
while in the remaining cases the 24/im positions were 
used as the target coordinates. An aperture correction 
factor appropriate f or a stellar photosphere was taken 
from Gordon ct al. (j2007[ ). After identifying stars with 
70/im excess, we recalculated their fluxes with an aper- 
ture correction appropriate for a 60 K blackbody (char- 
acteristic temperature of our disks). There are eleven 
images where bright nearby sources contaminated the 
photometry of our targets (see Notes in Table [5]) . In 
order to remove the contribution of these background 
objects, we fitted a PSF to them and then subtracted 
their emission (the PSF was constructe d based on the 
method described bv lGordon et al.ir2007[ ). For stars sur- 
rounded by bright extended nebulosity in the 70/im im- 
ages (HD 36248, HD 143840, HD 185053, HD 218980), no 
photometry was derived. Due to the resampling of 70/im 
mosaics, the noise between the adjacent pixels became 
correlated. In order to take into account this effect, in the 
course of internal photometric error estimation, we fo l- 
lowed the method described bv iCarpenter et all ()2008D . 
The resulting 70/im photometry is presented in Table [2) 

Positional offsets between the obtained centroids of the 
detected sources and their 2MASS positions as a function 
of signal-to-noise ratio {Fjo/ct^'q) measured at 70/im are 
displayed in Figure[31 Most of our detected 70/im sources 
are within 1.7" (^Icr pointing reconstruction accuracy at 
70/im, MIPS Instrument Handbook) of the correspond- 
ing stellar positions. It is worth noting that even in the 
case of HD 14691 - which shows the largest angular offset 
- the measured flux density is in good agreement with the 
predicted photospheric flux density at 70/im, conflrming 
the association between the 70/im source and the star. 
It is thus probable that all of the detected sources are 
likely to be associated with our target stars. 

Figure S] shows the histogram of the signiflcances of 
the differences between the measured and predicted pho- 
tospheric flux densities, deflned as (F70 — Pro)/cr7o', for 
all of our 78 targets (nebulous objects were excluded). 
The peak around zero signiflcance level corresponds to 
either stars that have not been detected at this wave- 
length or objects that show pure photospheric emission. 
A Gaussian fit to this peak yields a mean of -0.09 and 
a dispersion of 0.87, and these parameters are in good 
agreement with the expectations (mean of 0.0 with dis- 
persion of 1.0). 

The total noise in the photometry was derived as the 
quadratic sum of the internal and the absolute calibra- 
tion uncertainties (7%, see the MIPS Data Handbook). 

MIPS 160.— MIPS 160/im observations suffer from a 
spectral leak resulting in a ghost image at a certain 
offset from the nomin al position of the original target 
(e.g. IStansberrv et al.l |2(307|). According to the MIPS 
Data Handbook (ver. 3.2) the ghost is bright enough 
to appear above the confusion noise only for sources 
brighter than 5.5 mag in J band. For those of our tar- 
gets (10 objects, see Table[2]for details) that are brighter 



than this lirn i t, we applied the procedure proposed by 
iTanner et all (|2009l ) to minimize the contamination of 
the leak image. We downloaded from the Spitzer Archive 
24/im and 160/im BCD level data fo r 13 bright stars 
listed in Table 2 of ' Tanner eFall (2009) that have no ex- 
cess in MIPS bands. These data were processed iden- 
tically to our F-stars observations (see Sect. I3.1.ip . On 
the flnal images we determined the average offset of the 
ghost images from the position of the stars (obtained 
at 24/im images), yielding 5 and 0.2 pixels in the X 
and Y direct i ons, r espectively, in good agreement with 
ITanner et all (j2009D . Then we compiled the typical PSF 
of the leak image using the ghost images of the bright- 
est four stars. As a flnal step, utilizing the leak PSF 
and the estimated offset position of the ghost image, we 
applied an iterative cleaning method to remove the leak 
from the 160/im images of our brightest F stars. We 
performed aperture photometry to obtain flux measure- 
ments on MIPS 160/im images. The aperture was cen- 
tered on the source positions obtained on the 24/im im- 
ages. The aperture radius was set to 32", the background 
level was estimated using an iterative clipping procedure 
in an annulus extending from 64" to 128". We used an 
aperture correction f actor corresponding to a 50 K tem- 
perature blackbody (jStansberrv et al.ll2007D . For those 
stars where extended nebulosity or a nearby bright source 
contaminated the area of the aperture we give no pho- 
tometry at 160/im in Table [21 In some cases, nearby 
bright sources appear in the background annulus. To 
minimize their influence, we placed a 4 pixel radius mask 
on these sources before we estimated the sky level (see 
Table [2] for the affected sources). 

MIPS SED— Two of our targets, HD 50571 and 
HD 170773 were observed in MIPS SED mode as well. 
The spectra were extracted from the sky-subtracted mo- 
saic images using a flvc-pixel-wide aperture. Aperture 
correction was needed to correct for flux losses. Although 
HD 170773 is shghtly extended at 70/im (see Sect. 113)) . 
our studies based on smoot hed Spitzer TinyTim mod- 
els of the MIPS SED PSF (|Kristl 120021 iLu et al.l [2001 
showed that the difference in aperture corrections be- 
tween a point source and a point source convolved by a 
Gaussian with FWHM of 10" is insignificant. Thus we 
applied apertur e correction factors valid for point sources 
(jLu et al.ll2008D for both stars. We discarded the longest 
wavelength part of both spectra (A > 90/im) because of 
their very low signal to noise. As a final step of the data 
processing, we scaled our spectra to the photometry ob- 
tained in the 70/im band as follows (the MIPS SED and 
MIPS photometric observations were performed on the 
same day). We extrapolated the spectra towards both 
shorter and longer wavelengths based on IRS and 160/im 
photometry, in order to completely cover the transmis- 
sion curve of the MIPS 70/im filter. We computed syn- 
thetic 70/im photometry and derived the ratio of the 
synthetic to the measured flux density. According to 
the resulting ratios, the MIPS SED of HD 15745 and 
HD 170773 were scaled down by 1.11 and 1.16, respec- 
tively. 

3.2. IRS 

The reduction of the IRS observations started with in- 
termediate droopres products of the SSC pipeline (ver- 
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sion S15.3). We u sed the SMART reduction package 
(|Higdon et al.|[200^ . in combination with IDL routines 
developed for the FEPS Spitzer Science Legacy program 
([Meyer et al.l [2006) . We first subtracted the pairs of 
imaged spectra acquired along the spatial direction of 
the slit to correct for background emission, stray light, 
and pixels with anomalous dark current. We replaced 
bad pixels by interpolating over neighboring good pixels. 
Spectra were extracted from the background-subtracted 
pixel-corrected images using a 6-pixel fixed-width aper- 
ture in the spatial direction centered at the position de- 
rived by the FEPS Legacy program. After extracting 
the spectra for each order, nod, and cycle, we computed 
a mean spectrum for each order and as uncertainty we 
quoted the Icr standard deviation of the distribution of 
the fiux densities measured at the given wavelength. We 
converted our output spectra to flux units by apply- 
ing the spectr al response function d erived by the FEPS 
Legacy team ([Bouwman et al.ir2008D . and propagate the 
calibration error into the quoted uncertainties. 

Using the filter transmission curve of the MIPS 24/im 
band we computed synthetic photometry from the IRS 
spectra {IRS2A) for all of the targets. Figure [5] shows 
a comparison between the synthetic IRS 24/im and the 
measured MIPS 24/zm flux densities (^^24)- The calibra- 
tion of the two instruments shows good agreement in 
general, the mean of the IRS24/ F24 ratios is 1.02 with 
a dispersion of 0.06. The largest discrepancy between 
the two instruments (^30%) was found at HD 185053, a 
source surrounded with extended nebulosity. 

3.3. Additional data 

For the targets where the Spitzer observations pointed 
to the existence of excess emission (Sect. 14. 1|) we col- 
lected additional infrared and submillimeter data from 
the literature. IRAS 60/im and 100/im flux densi- 
ties and their uncertaint ies were ta ken from the IRAS 
FSC (M oshir et al.|[l989l ). [Williams fc Andrews! (( 20061 ) 
detected three of our targets (HD 15115. HD 127821, 
HD 206893) at 850Min usin g the JCMT/SCUBA instru- 
ment, while [Nilsson et aD ([2010) observed three sample 
stars (HD 17390, HD 30447, HD 170773) at 870fim with 
the LABOCA/APEX instrument. 

Observations obtained with ISOPHOT, the photome- 
ter on-board the Infrared Space Observatory were avail- 
able for ten sources. We observed three northern stars 
with the IRAM 30-m telescope at millimeter wave- 
lengths. Moreover, for some targets we performed optical 
spectroscopy. In the following we describe the details of 
the data processing. 

3.3.1. ISO/ISOPHOT 

Ten out of the 27 objects with confirmed infrared ex- 
cess were observed with ISOPHOT as well. For eight 
of these ten st a rs, th ere are published ISOPHOT fluxes 
in I Moor et al.l ([2006D . For the remaining two objects, 
HD 151044, HD 213617, we processed the ISOPHOT 
data in the same way as described in that paper and 
the result s are given i n Tab le [Sj HD 170773 was in- 
cluded in I Moor et ahl ([2006D . however, the photomet- 
ric results were extracted assuming a point-like source. 
The analysis of our MIPS data showed that the disk 
around HD 170773 is spatially extended with an extent 



of '-^lO" (Sect. 14. 3p . Assuming that the spatial ex- 
tent of the disk is similar at 60/im and 90/im we re- 
analyzed the JSOPffOT data by convolving the appropri- 
ate ISOPHOT beam profiles by a Gaussian with FWHM 
of ^10" and then using this new profile in the course 
of the flux extraction. The new photometry is given in 
Table [3 

3.3.2. IRAM/MAMB02 

We observed three stars (HD 15745, HD 25570, 
HD 113337) at 1.2 mm using the IRAM 30-m telescope at 
Pico Veleta with the 117-element MAMB02 bolometer 
(proposal ID: 195/07, PI: A. Moor). None of them have 
been observed at millimeter wavelengths before. The ob- 
servations were carried out between December 2007 and 
April 2008, using the standard on-off observing mode. 
The targets were always positioned on pixel 20 (the most 
sensitive bolometer pixel). Observations of Mars were 
used to establish the absolute fiux calibration. The value 
of the zenith opacity at 1.2 mm ranged between 0.17 and 
0.48 during our observations. We performed 3, 3 and 11 
of 20 minutes long ON-OFF scan blocks for HD 15745, 
HD 25570, HD 113337, respectively. Wc utilized the 
MOPSIC software package (R. Zylka) to perform the 
data processing using the standard scripts developed for 
the reduction of on-off observation data with sky noise 
subtraction. None of our targets were detected at 3ct 
level. The obtained flux densities and their uncertainties 
are quoted in Table [3| 

3.3.3. Ground based spectroscopy 

In July and August 2009 we obtained new high- 
resolution optical spectroscopy for stars in Table [H us- 
ing the 2.3-m telescope and the Echelle spectrograph of 
the Australian National University. The total integration 
time per object ranged from 30 s to 1800 s, depending on 
the target brightness. The spectra covered the whole vi- 
sual range in 27 echelle orders between 3900 A and 6720 
A , with only small gaps between the three reddest or- 
ders. The nominal spectral resolution is A/AA w 23 000 
at the Ha line, with typical signal-to- noise ratios of about 
100. 

All data were reduced with standard IRAlQ tasks, 
including bias and fiat-field corrections, cosmic ray re- 
moval, extraction of the 27 individual orders of the 
echelle spectra, wavelength calibration, and continuum 
normalization. ThAr spectral lamp exposures were reg- 
ularly taken before and after every object spectrum to 
monitor the wavelength shifts of the spectra on the 
CCD. We also obtained spectra for the telluric standard 
HD 177724 and lAU radial velocity (RV) standards /3 Vir 
(sp. type F9V) and HD 223311 (K4IH). 

The spectroscopic data analysis consisted of two main 
steps. First, we measured radial velocities by cross- 
correlating the target spectra (using the IRAF task fxcor) 
with that of the RV standard that matched the spectral 
type of the target - /? Vir was used for the F-type targets, 
HD 223311 for the lone late G-type target (SAO 232842). 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Each spectral order was treated separately and the re- 
sulting velocities and the estimated uncertainties were 
calculated as the means and the standard deviations of 
the velocities from the individual orders. For most of 
the targets, the two lAU standards gave velocities within 
0.1-0.5 km/s, which is an independent measure of the ab- 
solute uncertainties. The equivalent width of the 6708A 
Li were measured with the IRAF task splot. The pro- 
jected rotational velocity of t he targets was dete rmined 
via fitting theoretical models (jMunari et al.l[2005l ) to the 
observed spectra with the method. Table |4] summa- 
rizes the derived properties of the observed stars. 

4. RESULTS 

4.1. Identification of stars with infrared excess 

We used our Spitzer data to identify stars exhibiting 
excess at infrared wavelengths. First, the predicted pho- 
tospheric fiux densities were determined at the relevant 
wavelengths using the best-fit Kurucz models of the stars 
(see Sect. 12.11) . The average accuracy of the predicted 
far-infrared fluxes is estimated to be around 3%. The 
predicted flux densities of the stars for the MIPS bands 
{P24, P160) are listed in Tabled The significance 
level of the infrared excess was calculated in each photo- 
metric band using the following formula: 



where is the measured fiux density, P^, is the predicted 
stellar flux, while cr*°* is the quadratic sum of the uncer- 
tainty of the measured flux density and the uncertainty 
of the predicted flux density in the specific band. When 
Xiy was greater than 3 in any of the MIPS bands, the 
object was selected as a star with excess emission. Ap- 
plying this criterion, we identified 28 stars that exhibit 
IR excess. One of our targets, HD 199391, shows excess 
only at 160/im. Since at this wavelength the position 
measured on the 24/im image was adopted in the course 
of photometry, we cannot exclude the possibility that the 
excess emission is related to a nearby background source. 
Thus, HD 199391 was excluded from the list of stars with 
excess emission. Future observations with higher spatial 
resolution in far-infrared bands should reveal the true 
nature of the excess emission observed towards this ob- 
ject. Among the 27 remaining systems 15 show excess 
at 24/^m, all 27 exhibit excess at 70/im and we found ex- 
cess at 160/im in 17 cases. All of our targets that exhibit 
excess at 24/im also show excess at 70/xm as well. 

Although in most cases the shape of the obtained IRS 
spectra could be well fitted by the photospheric model of 
the specific target, 32 among the 82 spectra showed sig- 
nificant deviations. Out of the 27 objects where excess 
was indicated by the MIPS photometry, 25 also showed 
excess in the IRS spectra. In the case of HD 15060 
and HD 213429 the IRS spectra were consistent with the 
predicted photospheric emission and these two objects 
showed excess only at 70/im. Seven additional excess 
sources were revealed by the IRS data. In three of the 
seven cases (HD 143840, HD 185053 and HD 218980) the 
MIPS images show bright nebulosity around the stars 
with a spatial extent of 50-60" (see also Sect. I3.1.2p . 
At the distance of these stars the estimated angular ex- 
tents correspond to a size of 4000-7000 AU. Thus in these 



cases the observed emission is likely to be of interstel- 
lar, rather than circumstellar origin. The remaining four 
stars (HD 34739, HD 38905, HD 145371, HD 184169) are 
located close to background sources that are brighter at 
mid- and far-IR wavelengths than the original target. In 
these cases we investigated - using the beam profiles of 
the different IRS modules and the known position of the 
background sources (derived on the MIPS images) - the 
possibility that the observed excess emission was associ- 
ated with these bright nearby objects. We found that in 
all four cases the apparent excess is likely to be related to 
the nearby sources (Fig. |6] demonstrates the main steps 
of our analysis for HD 38905). As a consequence, these 
additional seven candidates were discarded from the fur- 
ther analysis. 

A significant fraction of our candidates that were se- 
lected based on previous IR observations turned out to 
be misidentifications in the light of the new Spitzer data. 
Several earlier works showed that IRAS-based debris can- 
didate lists are strongly contaminated by false identifi- 
cations (bogus debris disks) especially due to the low 
spatial resolution of IRAS observations at far-IR wave- 
lengths ((Kalas e t al. 2002; Moor et al. 2006; Rhce ct all 
l2007h . The most common reasons of the misidentifica- 
tion are: 1) confusion with background sources, where 
the IR emission is associated with a nearby object; 2) 
the presence of extended nebulosity where the IR emis- 
sion is of interstellar rather than of circumstellar origin; 
3) erroneous infrared photometric measurements. The 
last reason played an especially important role in the 
case of warm disk candidates. 

In total 27 stars have been identified where the ob- 
served excess emission may originate from circumstellar 
dust grains (see Table [6]). Nine out of the 27 stars ex- 
hibiting IR excess are new discoveries (see in Table |6]). 

4.2. Spectral features in the IRS spectra 

No prominent features have been identified in the IRS 
spectra of the 27 stars with excess emission. This finding 
is consistent with the results of previous IRS observations 
related to debris disks around solar-like and more mas- 
sive A-type stars: although many debris systems show 
significant excess in the wavelength range between 5 and 
35/im the maj ority of these systems do not possess spec- 
tral features (iChen et al.ll2006t ICarpenter et a"Ll l2009al: 
iLawler et afll2009l : [Morales et al.ll2009| l 

4.3. Spatially extended sources at 70^m 

Observations with the MIPS detector at 70/im revealed 
several (marginally) resolved debris disks even at rela- 
tively large distances from the Sun (e.g. iBrvden et al.l 
_2006; Su ct al. 2009). In order to identify sources with 
extended emission, in Figure [7] we plotted the ratio of 
the flux density measured in apertures with radius of 
8" and 18" as the function of the SNR obtained in the 
smaller aperture. Similar ratios derived fo r some known 
resolv ed debris disks (based on the list of IBrvden et al.l 
[2006h were also plotted. MIPS 70//m data for these ob- 
jects were downloaded from the Spitzer Archive and pro- 
cessed using the same method applied to our targets. As 
expected, objects with known extended disks show signif- 
icantly la rger flux ra tios. Utilizing the Spitzer TinyTim 
software (Krist"2002) and applying the smoothing proce- 
dure proposed bv , Gordon et al.. (2007.') we calculated the 
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PSF of a 60 K blackbody source that was used to derive 
the expected flux ratio plotted in Fig. [T) By comparing 
the obtained flux ratios to the expected one, we found 
that two of our sources - HD 50571, HD 170773 - show 
significant (>3cr level) deviation and may be spatially ex- 
tended at 70/xm. In order to model the observed profiles 
we convolved the PSF with different Gaussian profiles. 
The minor- and major-axis FWHMs as well as the posi- 
tion angle of the major axis of the Gaussians were varied 
to find the best-fitting model profile. We found that the 
measured profile for HD 50571 can be well fitted using a 
Gaussian broadened by 9.5" only in one direction along 
a position angle of 91° (measured from north to east). 
In the case of HD 170773, a convolution by a Gaussian 
with minor- and major-axis FWHMs of 9" and 10" with 
position angle of 110° provided the best fit. MIPS 70/im 
photometry for these sources, which takes into account 
the spatial extent, is given in Table [5] 

It is a common expectation that the rotational spin 
axis of the star is aligned with the orbital spin axis of the 
planetesimals. HD 50571 has an u nusually high projecte d 
rotational velocity of 60kms~^ (jHolmberg et al.ll2007D . 
When comparing this value to the vsini of 1566 stars 
with similar effective temperat ures (Teg = 6480 ± 10 K) 
included in the same catalog ()Holmberg et al.ll2007D we 
found that the projected velocity of HD 50571 is higher 
than 99.4% of that of the catalog sample. This implies 
a very high inclination of the spin axis and consistent 
with our finding that the disk around this star may be 
seen nearly edge-on (it is resolved only in one direction). 
The derived size of the debris disk around HD 170773 
indicates that it might be seen close to pole-on. It is 
interesting that the projected rotational velocity is re- 
markably high (see Table 0]) for this orientation. 

4.4. Modeling the observed infrared excess 

The infrared excess emission of F-type main-sequence 
stars is generally believed to be attributed to the opti- 
cally thin thermal emission of second generation circum- 
stellar dust grains heated by the central star. Since no 
resolved images are available for our targets, the funda- 
mental parameters of these dust disks have to be derived 
by modeling the SEDs of the systems. Four newly dis- 
covere d disks have alrea dy been modeled in a previous 
paper (|Mo6r "etal|[2009t ). For the remaining 27 stars ex- 
hibiting excess emission we compiled the SED from the 
data listed in Sect. [3] For the fitting process the IRS 
spectra were sampled in 11 adjacent bins (for the center 
and width of the bins, as well as for the obtained flux 
densities see Table [S]) . We used a simple model to char- 
acterize the disk properties based on the excess emission. 
We assume that the dust grains are distributed around a 
single radius and we adopt the same temperature for all 
particles within this ring. Then the excesses are fitted 
by a single temperature modified blackbody, where the 
emissivity is equal to 1 at A < Ao and vary as (A/Ao)~^ 
at A > Ao wavelengths. We used this modified blackbody 
model in order to account for the falloff in the emission 
spectrum at longer wavelengths, that is faster than in 
the ca se of a blackbody. Following IWilliams &: Andrewsl 
(|2006f) we fixed Aq to 100/xm consistent with the lack of 
spectral features in the IRS spectra which indicates a 
relatively large average grain size. Since /3 cannot be 
reliably determined in all systems due to the lack of 



long wavelength data, we decided to estimate a char- 
acteristic /? value based on those disks where the excess 
emission at A > 50/.tm was measured with > 3 in 
at least 4 different bands and were successfully detected 
at A > 100/im as well. A Levenberg-Marquardt algo- 
rithm was used to fit the model to the measured data 
and an iterative method was used to compute and apply 
color corrections for t he photometric da ta during the fit- 
ting process (see e.g. iMoor et al.ll2b06D . Figure |8] shows 
the obtained /? parameters for the selected 11 disks. A 
characteristic (3 value of 0.7, derived by computing the 
weighted average of these values, seems to represent very 
well the whole sample. Thus, in the following we fixed 
(3 equal to 0.7 and Aq to lOOyitm and repeated the fitting 
process for all debris systems. In two cases when the 
excess was detected only in the MIPS 70^m band, we 
determined the highest possible dust temperature which 
was still consistent with the IRS data. The derived dust 
temperatures (Tdust) and the reduced chi-square values 
of the best fits are quoted in Table [6] We note that us- 
ing a simple blackbody model instead of the modified 
blackbody would change the dust temperatures less than 
their formal uncertainties, but would yield worse fitting 
for the long wavelength data in most cases when we have 
A >100/xm measurements. 

Most of our SEDs can be fitted well using the sim- 
ple model described above. However, in several cases 
systematic deviations between the model and measured 
fluxes can be seen, especially at shorter wavelengths, 
where the model typically underestimates the observed 
excess. For example, for HD 16743 and HD 192758 the 
high reduced chi-square values indicate poor model fit. 
It is a general trend among the deviating cases that fit- 
ting only the IRS spectra yields higher dust tempera- 
ture estimates than fits based on the long wavelength 
photometric points. We also found that the extrapola- 
tion of IRS-based fits to longer wavelengths underpre- 
dicts the excess measured in the IRAS, ISOPHOT and 
MIPS bands. Several authors reported similar findings 
in different debris disk samples ([H illenbrand ct al. 200q; 
iCarpenter et al.ll2009"at iMoraleset al. 2009). One possi- 
ble explanation for the observed discrepancy is that dust 
grains in these systems are distributed in two spatially 
separated rings similarly to our Solar System, where the 
majority of dust grains are thought to be co-located with 
the main asteroid and with the Kuiper belt. Assuming 
that the dust is concentrated in two distinct narrow rings, 
we used a two-component model, where grains in the 
warmer component act like blackbodies, while the emis- 
sion of the outer ring can be described by the modified 
blackbody as defined above. To decide whether a single 
or a two-component model should be used for a certain 
target, we used a variant of the Akaike Information Crite- 
rion, t he so-called AICu, proposed bv lMcQuarrie fc Tsail 
(IT998I) . The value of AICu can be calculated as: 



AICu = In 



SSEk 2(fc + l) 



fc-2' 



(2) 



where n is the number of observations, k is the number 
of parameters in the model, while SSEk is the usual sum 
of squared errors. Besides the fact that the Akaike Infor- 
mation Criterion take into account the goodness of fit, it 
penalizes the usage of unnecessary additional model pa- 
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rameters. This test can be used to rank the competing 
models, from which the best one gives the lowest AICu 
value. We found the two- component model to be better 
in five cases (HD 15115, HD 15745, HD 16743, HD 30447, 
HD 192758). The derived parameters of these five disks 
are presented in Table [T] The SED and the best-fit mod- 
els for each of our targets are plotted in Figure M 

The fractional luminosity of the disks was computed 
as /dust = idust/iboi- The integrated dust emission was 
derived based on the fitted model while the star's lumi- 
nosity was calculated from the best-fit Kurucz model. 
We estimated the radi us of the dust ring (or ring s) using 
the following formula (jBackman fc Parescelll993[ ): 



Rdv 



AU 



Lr. 



/'278K^ ^ 

V ^duE 



(3) 



Because this formula assumes blackbody-like grains, the 
resulting i?dust corresponds to a minimum possible ra- 
dius. The obtained fundamental disk properties are 
listed in Tables [6] and [T] 

We made several simplifying assumptions in the ap- 
plied model. Collisions among planetesimals in debris 
disks produce a coUisional cascade, in which collisions 
gradually grind large bodies into smaller ones that are 
removed by radiation forces. This process is thought to 
result in a characteristic dust grain size distribution and 
copious amounts of small dust (Wyatt 2008, and refer- 
ences therein) . Even if planetesimals are distributed in a 
narrow ring, the radiation pressure pushes the smallest 
grains into more eccentric orbits extending the dust disk 
outward. Moreover, the planetesimal ring(s) can be ex- 
tended. Our model assumes relatively large grains that 
act like a blackbody at A < 100/im located in one (or 
two, see above) ring(s) of dust. Both the existence of 
smaller dust grains - that are ineffective emitters and 
therefore have higher temperature than large grains at 
the same radial distant from the star - and the finite 
radial extension lead to multi-temperature distribution. 

Two of the multiple temperature disks from Table [7] 
were resolved or marginally resolved in scattered light. 
Note, however, that different observation techniques are 
sensitive to different populations of dust grains. The 
scattered light images are expected to trace very small 
grains which may show different spatial distribution than 
those grains which dominate the mid- and far-IR emis- 
sion of the disk (Wyatt 2006), because a significant 
fraction of small grains could be blown out from the 
system outside t he planetesimal rin g. U sing corono- 
xaph ic images, iKalas et"al] ()2007aD and iDebes et al.l 
20081) successfully resolved a very extended circumstel- 
lar disk around HD 15115 at optical and near-IR wave- 
lengths, revealing a strongly asymmetric disk structure. 
The lobes of the disk can be traced inward to ^^31 AU. It 
could be consistent both with our single/two- ring mod- 
els where the outer ring is located at ^^40 AU. HD 15745 
was also resolved in scattered light using the Advanced 
Camera fo r Surveys aboard the Hubble Space Telescope 
(|Kalas et al . 2007b). The circumstellar disk is detected 
between ~128-480AU radius. The detection at the in- 
ner part of the disk is limited by PSF subtraction ar- 
tifacts, thus this image does not provide further infor- 
mation about the disk morphology in the inner regions. 
Our model with an outer ring at 21 AU disk is not incon- 



sistent with this result. Since the existing data do not 
allow to localize the warm dust unambiguously, in the 
following analysis of the five disks in Table [7] we assume 
that the warm emission originates from an inner ring. 

Two additional disks in our sample were marginally 
resolved in the 70/im images (Sect. 14. 3|) . The measured 
broadening suggests a dust structure size of 320 AU and 
370 AU for HD 50571 and HD 170773, respectively. Al- 
though both stars harbor relatively cold and extended 
dust rings, the minimum diameter of 132 AU and 146 AU, 
derived from our simple model, are significantly lower 
than the measured extensions. One possible explana- 
tion would be a dust structure similar to that of the 
F5/F6 V star HD 1813 2 7 (no t included in our sample), 
where iSc hneider et al.l (|2OO60 discovered a dust ring us- 
ing NICMOS coronagraphic observations. The extension 
of this ring (^^86 AU) significantly exceeds the modeled 
dust radius of 22 AU derived fro m the SED of t he ob- 
ject assuming blackbody grains (jSchneider et al.l [2006) . 
They propose that a large amount of small dust particles 
in the ring can explain the observed discrepancy since 
small grains are hotter than large grains at the same 
location. An analogy with Vega can provide another ex- 
planati on. Using M IPS observations that resolved the 
source. ISu et al.l (|2005i ) found that the radius of the disk 
at 70/im exceeds significantly the size of the disk seen at 
submillimeter wavelengths. The discrepancy probably 
comes from an extended cloud of small particles blown 
away from a planetesimal ring by the radiation pressure 
of the star. Due to the marginal resolution of the Spitzer 
70/im observations of HD 50571 and HD 170773, we can- 
not decide between the two scenarios. In the further 
analysis of these two objects we adopt the parameters 
fisted in Table El 

Using the derived disk parameters in Table [6] and [7] we 
estimated both the Poynting- Robertson (tpr) and the 
collision timescales (tcou) for grains with radii ranging 
between the grain size corresponding to the blowout limit 
and 1000/im. The blo wout lim it was computed using the 
equation presented bv iHillen brand et al. (2008): 



0.52 



2.5 gem L^/Lq 
p (T,/5780)' 



(4) 



where p is the density of the grain (assumed to be 
2.7gcm~^). The Poynting-Robertson timescales were 
estimated based on equation 14 in Backman & Parcs(3 
(|1993| ) , while the collisional timescales were computed us- 
in g: the semi-empirical formulae (equations 7-8) derived 
bv iThebault fc Augereaul ()2007D . We found that for our 
disks Tcou <C TpR for all grain sizes and the obtained 
timescales are short with respect to the ages of the stars, 
implying that the grains have second generation ('de- 
bris') nature and their evolution is mainly governed by 
collisions. 

In a debris disk the mutual collisions continually grind 
down the larger planetesimals into smaller fragments 
that can be removed by the Poynting-Robertson drag 
and by the radiation pressure. Since the estimated colli- 
sional timescale in our disks is significantly shorter than 
the PR-drag lifetime, the dust removal processes in these 
systems may be collisionally dominated. In such a disk, 
the frequent collisions shatter dust grains more rapidly to 
sizes below the blowout limit before the effect of the PR 
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drag can be manifested. Due to the stellar radiation pres- 
sure, grains smaller than the blowout limit are ejected 
on a short timescale, while somewhat larger grains are 
pushed into a more eccentric orbit. Thus, in a coUi- 
sionally dominated disk, the dust grains extend outward 
from their birth ring, where the planetesimals are lo- 
cated. However, since the lifetime of blowout grains is 
significantly shorter than the normal grains' lifetime, it 
is a plausible assumption that the dominant part of the 
dust mass is co-located with the parent planetesimals. 
Thus in the further analysis we assume that the derived 
radii of the dust rings can be considered as the size of 
the underlying planetesimal belts as well. 

4.5. Age determination 

In order to estimate the age of the 27 debris systems, 
we use the the following general strategy. The most 
accurate and reliable dating can be derived via clus- 
ter membership. Thus, if a specific target can be as- 
signed to a stellar kinematic group then we adopt the 
age of the group for the star. Beside the nine previ- 
ously known cluster members, we classify five new young 
moving group members among our targets (for more de- 
tails related to the new assignments see below). For 
field stars, isochrone fitting combined with diagnostics 
of rotation-driven activity indicators and lithium content 
(especially that of late-type companions) are considered 
in the age estimates. In the analysis of chromospheric 
an d coronal activity indicators, th e calibration derived 
by iMamajek Sz Hillenbrandl (|2008[ ) is applied whenever 
it was applicable (i.e. the B — V color indices fall in the 
appropriate range). In the case of lithium content, our 
age estimates are based on comparisons with the dis- 
tribu tion of similar propertie s in well-dated open clus- 
ters (iSestito fc Rand ich'^OO l) and y oung moving groups 
(|Mentuch et al.ll200 8: da Sil va et al.ir2009l ). Pre-main se- 
quence evolutionary models are used if the lithium con- 
tent or the activity indicators measured in the target or 
in its late- type companion indicated that the specific sys- 
tem may be in a pre-main sequence evolutionary stage. 
We use literature data for several stars. Among the 27 
debris systems, 13 are younger than 100 Myr. Table |8] 
summarizes the age-related data for stars with debris 
disks. 

HD 3670. — The derived galactic velocity components of 
HD3670, U=-12.6, V=-22.6, W=-4.6kms-\ are con- 
sistent with the character istic motion of the 30 Myr old 
Columba association fsee iTorres et al.ll2008[ ). IID3670 
has a ROSAT counterpart with fractional X-ray lumi- 
nosity of log(Lx/Lboi) — —4.39, which is comparable with 
the fractional X-ray luminosity of those stars with similar 
spectral type in the Columba association and supports 
the youthfulness of the star. 

HD 15745. — Based on our new radial velocity measure- 
ments (Table El, we recomputed the galactic space veloc- 
ities of this star, obtaining —10.4, —15.3, — 7.9kms~^ for 
the U, V, W components, respectively. This space veloc- 
ity corresponds well to the c haracteristic space m otion of 
the /3 Pic moving group (see iTorres et ahfeOOSD . 

In the framework of a recent high-resolution spec- 
troscopic survey focusing on optical counterparts of X- 
ray sources. iGuillout et al.. (,2009. ) discovered several new 



post-T Tauri stars located on the northern hemisphere. 
One of them, BD-t-45° 598 - which is offset by ~9° from 
HD 15745 - shows very similar proper motion and ra- 
dial velocity {iiacosS = 44.7 ± 1.1 mas, us = —44.3 ± 
1.0 mas, = —0.77 ± 0.97kms^^) to that of our tar- 
get (/iQcos(5 = 45.8 ± 0.6 mas, = —47.9 ± 0.5 mas, 
V, = +2.5 ± 3.3 k ms"^). BD-h45° 598 was classified as 
a Kl type star bv^Gui llout et all (|2009| ). The measured 
lithium equivalent width and the fractional X-ray lumi- 
nosity of this star are consistent with the similar prop- 
erties of the known /3Pic members (see Fig. [10] a, b). 
Assuming that BD-|-45° 598 also belongs to the (5 Pic 
moving group, a kinematic distance of 70 pc can be esti- 
mated for it, obtaining —10.4, —16.2, — 8.1kms^^ as its 
galactic space motion. The position of the two stars on 
the color-magnitude diagram of the /? Pic moving group 
also confirms their membership status (see Fig. [10] c). 
The fact that HD 15745 harbors a debris disk with very 
large fractional dust luminosity also supports this assign- 
ment. 

Therefore, we propose that both HD 15745 and 
BD+45° 598 are new members of the /? Pic moving group 
and we adopt an age of 12 Myr for HD 15745. 

HD I6743 — Both the proper motion [iiacosS — 73.12 ± 
0.27 mas, jis — 49.65 ± 0.3 mas) and the trigonomet- 
ric parallax (tt = 16.99 ± 0.31 mas) of HD 16743 are 
in good agreement with the corresponding astrometric 
properties of HD 16699AB {pL^cosS = 72.33 ± 0.86 mas, 
^is = 48.56 ± 0.87 mas, tt = 16.54 ± 0.99 mas), that is 
also a multiple system itself with a separation of 8'.' 7 
(HD16699-I-SAO 232842). The similarities of the mea- 
sured radial velocities of the three stars (see Table [J]) 
also confirm that they may form a wide multiple sys- 
tem. This offers a good opportunity to improve the age 
determination of HD 16743 by combining the result of 
different age diagnostic methods for the three members 
of the system. 

The ROSAT source J023845.4-525710 is located close 
to both HD 16699 (with separation of 12'.'6) and 
SAO 232842 (3"9). This X-ray source is also present 
in the XMM-Newton slew s urvey Source Cata logue 
(XMMSLl J023845. 1-525708, ISaxton et al.][2008h . lo- 
cated 1" away from the position of SAO 232842. Due 
to the better positional accuracy of the XMM catalog, 
the latter data make it clear that the X-ray source cor- 
responds to SAO 232842. The high fractional X-ray lu- 
minosity of the source, L^/L\jo\ = —3.32, is compara- 
ble with that of stars with similar spectral type in the 
Pleiades and in young nearby moving groups confirming 
a young age for this object. 

The age estimates based on the lithium abundance 
measured in HD 16699 and SAO 232842 are some- 
what controversial. The high lithium abundance in 
SAO 232842 suggests a very young age, it exceeds the 
upper envelope of the distribution of lithium equivalent 
width measured in Pleiades stars and consistent with the 
lithium content of late-type G stars belonging to very 
young (<40Myr old) stellar kinematic groups and open 
clusters. On the other hand, the lithium equivalent width 
measured in HD 16699 suggests an older age, since it is 
somewhat lower than for stars with comparable effective 
temperatures in the Hyades and Coma Berenice open 
clusters (age ~ 500 - 600 Myr). 
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Figure [TT] shows the positions of the three components 
in the H-R d iagram overplotted by isochrone models from 
iSiess et alj ()2000[ ) for ages between 10 Myr and 50 Myr 
for a metallicity of Z=0.014. Assuming, primarily based 
on the characteristics of SAO 232842, that the system is 
in a pre-main sequence evolutionary stage, the positions 
of the three components suggest an age between 10 Myr 
and 50 Myr. 

HD 16743 and HD 16699AB belong to a very wide mul- 
tiple system with a minimum separation of ~12700AU. 
Gravitational perturbations due to the subsequent en- 
counters with other stars and giant molecular clouds dur- 
ing the galactic orbit would detach such wide binaries. 
Thus the existence of this system is also an additional 
argument concerning the young age. 

Taking into account age estimates obtained from the 
different dating methods we propose an age of 10-50 Myr 
for the HD 16743 system. 

Both HD 16699 and SAO 232842 are detected and their 
fluxes are derived via PSF fitting on the MIPS 24/im 
image of HD 16743. The obtained flux densities of 
16.4±0.7mJy and 18.8±0.8mJy correspond well to the 
predicted photosperic fiux densities in this band for both 
targets. 

HD 24636. — The derived galactic space motion of 
HD 24636 (see Table |4]) is consistent with the character- 
istic space moti on of the 30 Myr ol d Tucana-Horologium 
association fsee lTorres et al.ll200 8^. HD 24636 is located 
quite close to two other Tucana-Horologium stars from 
our sample (HD 32195, HD 53842), and all of these ob- 
jects are within a sphere 14 pc across. Interestingly, 
all three stars exhibit relativel y warm infrared excess 
(Tdust > 90 K, see Table [6] and iMoor et al.l[2009( ). 

HD 25570.— \Egged (il982[ ) classified H P 25570 as a star 
that belongs to the Hyades group. iPerrvman et al.l 
(|i998) found that HD 25570 is only just outside their 
3tT membership criteria and speculated that may be an 
example of an object just moving away from the Hyades 
but still close to the tidal radius. For the further analysis 
we consider HD 25570 as a Hyades member and adopt an 
age of 625 Myr for the star. 

HD 36968. — Both the galactic space motion 
([U,V,W] = [-14.8,-6.6,-8.5]km-i) and the galac- 
tic position of HD 36968 ([X,Y,Z] = [-51,-108,-73] pc) 
are consistent with those of the recently discovered 
~ 20 Myr old O ctans Association (see Fig. 17 in 
[Torres et alll2008D . HD 36968 is the first member of 
this association where the existence of a debris disk 
was revealed by our results. The very high fractional 
luminosity of this debris disk is also consistent with this 
young age. 

HD 50571 — is an F7-type star th at shows both m oder- 
ate chromospheric (i?yj^=— 4.55, IGrav et aI]|2006D and 
coronal activity (log 5.27). Its galactic motion 

(U = -16.6±0.3, V=-22.3±0.7, W= -4.4±0.3kms-i) 
as well as its galactic position are consistent with those 
of the " B3" subgroup of t he Local Associatio n, implying 
an age of 300±120Myr (jAsiain et al.lll999D. Us ing an 
independent age estimating method. iRhee et al.l (|2007[ ) 
also derived 300 Myr as the age of the object. 



HD 113337 — has a late- type companion with a sep- 
aration of ~120". The companion's proper motion 
{fj-acosS = —174 ± 4mas, fig = 26 ± 1 mas) is consis- 
tent with that of HD 113337 {fiacosS — — 171.7±0.3nias, 
fis — 25.4 ± 0.2 mas). B ased on me dium resolution spec- 
troscopic observations, IReid et all (j2007i ) classified the 
companion as an M3.5 dwarf with strong emission 

(log ^^^=—3.31). The companion is also included in the 

2MASS and SDSS (Sloan Digital Sky Survey) catalogs. 
Utilizing the spectroscopic a nalysis of a l arge n umber of 
M-type stars from the SDSS. I West et al.l (|2008[ ) provided 
characteristic median SDSS/2MASS colors {r—i, i—z, z — 
J^J—H, H—K) for different spectral subclasses (see their 
Table 1). Comparing the measured colors of HD 113337B 
with the quoted median values, we classify i t as an M4- 
type star, confirming the result of lReid et al.l ()2007[ ). Us- 
ing the spectral class information we estimated the ef- 
fective tem perature of the star utilizing the formula of 
IGrav fc Co rballv (2009), T^s[K] = 3759 - 135a;, where 
X is the optically defined spectral type valid for optical 
types MO through L8 with x = for an MO type, x = 4 
for an M4 type star, etc. Adopting 3.5 for the value of 
X this equation yields ~3290K as the effective temper- 
ature of HD 113337B. The age of the companion is esti- 
mated by comparing its position in the H-R diagram with 
the evolutionary tracks compiled by Siess_et al. (2000) , 
yielding ^40±20Myr (see Fig. [T2|) . This age estimate is 
adopted for t he whole system ( a similar age estimate was 
proposed by IRhee et al.ll2007n . 

Aperture photometry for HD113337B on the 24/im 
image yields a flux density of 2.38±0.22mJy. After a 
conversion of the 24/im MIPS fiux to magnitude we can 
compare the measured Kg — [24] color of the star to the 
stellar photospheric colo rs determined from a study of 
nearby M-type stars by IGautier et al.l ()2007D . On the 
basis of this comparison (see Fig. fT^ HD 113337B shows 
>3ct excess at 24/im. The presence of circumstellar dust 
may explain the observed excess. However, since the sig- 
nificance of the excess detection is just above the 3cr level, 
confirmation at additional wavelengths with higher sensi- 
tivity and spatial resolution is desirable. These observa- 
tions can also help to exclude alternative explanations of 
the apparent excess, like the possibility of contamination 
by a background galaxy or the presence of an unresolved 
low mass companion. Assuming that HD 113337B hosts 
a debris disk, we note that the number o f known debris 
disks around M-type stars is very limited (iForbrich et al.l 
[200l and all sources in this small sample showing excess 
at 24/im are younger than 50 Myr. 

HD 205674- — Apart from the star's galactic space veloc- 
ity component toward the Galactic center, that differs 
from the average U velocity component of the AB Dor 
moving group more than 2tT, the other velocity and space 
components are in good agreement with the similar prop- 
erties of th is kinematic assem blage (see Table 1 and 
Fig. 21 in iTorres et al1[200l . HD 205674 fits weh to 
the locus of AB Dor stars in the color-magnitude dia- 
gram. It ha s an X -ray counterpart as well (see Table [5]) . 
IRhee et all (|2007[ ) proposed an age of 300 Myr. Since 
currently no reliable age-dating criteria for the youth of 
HD 205674 are available and the membership status is 
questionable because of the deviation of U space veloc- 
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ity component from the cluster center, we adopt an age 
range wh ich covers both the age of the moving gr oup (70- 
150 Myr, iLuhman et aLl2005fllbrres et aLl2008l) and the 
300 Myr. 

5. DISCUSSION 

Our investigation of 82 F-type stars with the Spitzer 
Space Telescope resulted in the detection of 27 debris 
disks, out of which 9 are new discoveries. In the following 
we analyze the parameters of these disks with special 
attention to disk evolution and host star properties. 

5.1. Metallicity 

Table [1] lists metallicity estimates for 24 out of the 
27 disk bearing stars. The average metallicity value in 
this sample is -0.09±0.09. For comparison we selected 
9138 stars from th e Geneva- Copenhagen Survey of So- 
lar neighbourhood (jHolmberg et al.l[2007| ) with effective 
temperature falling in the range spanned by our sample. 
Their average metallicity, -0.11±0.22, is in good agree- 
ment with the result for our stars. Narrowing the com- 
parison sample to stars with similar age range would not 
change the conclusion. Thus, our sample is similar to 
stars located in the Solar vicinity in terms of metallicity. 
This is in accordance with findings that the incidence 
of deb ris disk s does not correlate with stell ar metallicity 
(jBeichman et HI l200l [Greaves et al.ll2006[) . The lack of 
any such correlation may suggest that the formation of 
planetesimals is not sensitive to the metallicity in the 
protoplanetary disks. 

5.2. Multiplicity 

The effect of binarity on the pres ence of debris disks 
was studied bv iTrilling et all (|2007D . They found that 
the incidence of debris disks is ~50% in systems with 
small (<3 AU) or wide (>50 AU) separations, even higher 
than the corresponding value among single systems. In 
our 82-star sample there are 23 known multiple systems 
(13 have known separation). Three of the multiple sys- 
tems harbor debris disks. Two disks are associated with 
the widest binaries with separation >4400 AU, where the 
components practically can be regarded as isolated stars. 
It is interesting to note that in the case of HD 113337, 
the secondary component might also harbor a debris disk 
based on the 24^m image (Sect. |4?5|) . The third disk en- 
circles HD 213429 , which is a spe ctroscopic binary whose 
orbital solution (|Pourbaixl l2000f ) indicates the smallest 
known separation in our sample (1.8 AU), probably form- 
ing a circumbinary structure. Thus, all three binaries fall 
in groups of multiple systems where the incidence of de- 
bris disks is high according to ITrilling et al.l (j2007[ ). 

5.3. Disk temperature 

The dust temperature provided by our modeling 
(Sect. is a fundamental disk parameter, whose com- 
putation includes only a few assumptions and can be 
determined with confidence. In Fig. [14] we plot a his- 
togram of the derived dust temperature values. In those 
five cases when the SED was modeled by multiple dust 
rings (Table [7|) only the colder component was taken into 
account. Disks with an upper limit for the tempera- 
ture were also omitted. Most of the disks in Fig. [14] 
have temperatures falling in the range of 40-70 K (in 



the Solar System the Kuiper-belt exhibits similar tem- 
perature). A smaller sample shows temperatures 70 < 
Tdust <120K. Note that our complete sample of F-type 
stars contains also four disks of even higher temperature 
(135- 200 K) published in our previous paper ( Moor et all 
|2009() . While the age of the disk host stars belonging to 
the colder group ranges between 12 Myr and 3,100 Myr, 
all warmer disks encir cle stars with age <150 Myr. Us- 
ing the formulae 1-5 in lGrigorieva et all (|2007[ ). we com- 
puted the sublimation time for possible ice grains with 
radii range between the blowout limit and 10/im using 
the derived grain temperatures. Apart from the case 
of HD 24636 the timescale of sublimation is significantly 
longer than the age of the system. Thus the presence of 
icy grains cannot be excluded in most of our disks. Note, 
however, that in such disks the photo-desorption can re- 
move icy mantles o f grains more effectively than the sub- 
limation ([Gri gorieva et al.ll2007[) . Moreover, grains are 
thought to be produced in collisions that can provide 
enough heat to subliinate th e icy mantles of the particles 
(iCzechowski &: MannI 120071 ) . In HD 24636, the sublima- 
tion timescale is lower than the age of the system even 
for icy grains with a radius of a few centimeters. Taking 
into account the uncertainties of the timescale estimates 
and the derived parameters in this disk, the sublimation 
and the coUisional timescales for the smallest grains are 
in the same order. If the grains are icy in this system, 
this finding suggests that sublimation may have a role in 
the removal of small grains. 

5.4. Fractional luminosity 

Fractional luminosity is a fundamental observable pa- 
rameter of debris disks that is frequently used to char- 
acterize the amount of circumstellar matter. Fig. [15] 
shows the derived fractional luminosities for our disks as 
a function of the system's age. Although the dust frac- 
tional luminosity shows a large dispersion at any given 
age, a clear decline with time can be recognized. Follow- 
ing the evolution of a narrow debris ring co-located with 
a planetesimal belt in an analytical steady state coUi- 
sional evolution model, IWvatt et al.l (|2007al) found that 
the fractional luminosity of the disk varies with time as 
t^^. In the f ramework of a mor e realistic numerical 
evolution model, iLohne et al.l (|2008[ ) lifted some simplifi- 
cations of the above mentioned analytical model by tak- 
ing into account that planetesimal strength depends on 
their size and the complex grain size distribution close 
to the blowout limit. They predicted that the decay of 
dust lu minosity is proportional to t~°', with a = 0.3 
to 0.4. iKenvon fc BromlevI (|2008D dealt with the for- 
mation and evolution of an extended planetesimal disk 
using a numerical model. They found that the maxi- 
mum in the dust emission coincides with the formation 
of 1000-2000 km size planetesimals at the inner edge of 
the disk. Therefore, the fractional luminosity of a disk 
starts to rise when large planetesimals appear at the in- 
ner edge of the disk and stir the motion of smaller bod- 
ies, initiating their destructive collisions. After the rise 
to a peak brightness, a slow, roughly power-law decline 
is predicted, /dust oc t~" with a ~ 0.6 — 1.0. In Fig. US] 
we plotted the two extremes of these model predictions 
(a=0.3 and Q!=1.0). The distribution of the data points, 
in particular their upper envelope, seems to suggest a 
decay rate halfway between the two extremes. 
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The steady state evolution model of iWvatt et alj 
(|2007af) predicted that in the course of the evolution of 
a narrow debris ring, at any given age there is a maxi- 
mum fractional dust luminosity (/max), since originally 
more massive disks eat up their mass faster. Disks with 
fractional luminosity of 3> /max could be the result of 
a transient event that increases the dust product ion for 
a short period (for possible transie nt effects, see iWvatti 
120081) . Using their formula (eq 20. in lWvatt et al.ll2007aD 
and adopting their fixed model parameters (belt width: 
dr/r = 0.5, planetesimal strength: Ql, = 200 J kg^^, 
planetesimal eccentricity: e — 0.05, diameter of largest 
planetesimal in cascade: — 2000 km) we computed 
the /max values for our debris systems. By comparing 
the measured fractional luminosities (see Tables |6] and [71) 
with the calculated maxima we conclude that all our cold 
disks can be consistent with a steady-state evolutionary 
scenario within the uncertainties of the model. In reality, 
transient processes might be present but their influence 
is undetectable with this simple comparison. Note that 
this model uses a pre-stirred planetesimal belt (i.e. the 
stirring is init iated at f = ) and that in the numerical 
model of Loh ne et al.l ()2008[ ) the dust luminosity depends 
on the initial mass even at late evolutionary stages. Nev- 
ertheless, taking into account the delayed stirring or the 
dependence on initial disk mass do not affect our conclu- 
sions since both alterations would increase the /max ^t a 
specific age. 

Due to the decline in dust mass with time, high frac- 
tional luminosity debris disks are associated mainly with 
young stars. The rare exceptions tend to harbor very 
hot debris dust, in which the ejection of a small amount 
of transient dust can cause a significant rise in the frac- 
tional lumino sities (b ecause of the disk's proximity to 
the star). iMoor et al.l (jlQOQ) hypothesized that all stars 
with /dust > 5 • lO"** are younger than 100 Myr, there- 
fore a high /dust values can be used as an age indicator 
(see also in lZuckerman &: Sonsjl2004al but with a limit 
of /dust > 10"'^). Among our new discoveries there are 
two debris disks with fractional luminosity exceeding the 
limit of /dust > 5 • 10-^ HD3670 and HD 36968. Both 
objects are proposed to belong to young moving groups 
(age < 30 Myr) in agreement with the hypothesis. 

5.5. Disk radii 

In the course of disk radius estimates we assumed the 
dust grains to be confined to a narrow ring and that 
they interact with the stellar radiation as a blackbody. 
The resulting i?dust values correspond to minimum pos- 
sible radii and possibly underestimate the "true" radii, 
meaning that the real size of a specific debris disk could 
depart significantly from the derived -Rdust- However, if 
the disks are composed of similar dust grain populations 
then the differences between the real dust distribution 
and the assumed one would s hift the computed values 
in a similar way ()Wvattl[2008D . i.e. the relative radii of 
the disks are better constrained than the absolute val- 
ues. Thus, in the following analysis we assume that the 
derived values can be used to study general trends in the 
disk radii distribution. 

Destructive collisions between planetesimals can occur 
when the collision velocity exceeds a critical value that 
requires a dynamically excited (stirred) disk. In self- 



stirring models the formation of large planetesimals in 
collisional coagulation among smaller planetesimals nat- 
urally leads to the formation of a debris ring as well. 
These oligarchs can stir up the motion of the leftover 
smaller bodies initi alizing a collisional cascade. Accord- 
ing to the models of lKenvo"n~&: Br omlev (_2008) , the max- 
imum of the dust production via these collisions coin- 
cides roughly with the formation of ~1000km planetes- 
imals in the same region. Since the formation of such 
large bodies requires longer time at larger radial loca- 
tions, the site of the dust production in an extended 
planetesimal disk is thought to propagate outward. Sec- 
ular perturbations by giant planets - formed previously 
in the inner regions of the protoplanetary disk - can 
also initialize a co llisional cascade in a planetesimal disk. 
iMustih fc WvattI ([2009) concluded that planetary stir- 
ring can also eventuate in an outwardly propagating dust 
ring. In some regions the time-scale of this process can 
be even s horter than the growth time of ~1000km plan- 
etesimals (jMustin fc Wvattll2009[) . Stellar flybys can also 
initiate more energetic collisions in a planetesimal disk. 
However, such rare events are not likely to be responsible 
for large numbers of debris systems. 

Figure [TBI shows the derived radii of the dust rings as a 
function of age. The radii of the rings show large disper- 
sion at any given age. The data points seem to suggest 
an increase of the upper envelope of the distribution with 
increasing age. It is even more salient that while older 
systems (age > 100 Myr) harbor dust rings located at 
radii of >30AU, around younger systems there are sev- 
eral disks at radial location between 10 AU and 30 AU 
in a region where Saturn, Uranus, and Neptune orbit in 
the current configuration of our Solar System. The inset 
in Figure [T6l shows a comparison between the cumulative 
distribution of disk radii around stars with age <100 Myr 
and stars with age > 100 Myr (disks with lower limit for 
radius are not included) . This comparison suggests a sig- 
nificant difference between the two distributions. Based 
on a Wilcoxon test, the null hypothesis that the two sam- 
ples come from identical populations can be rejected on 
a 99.9% confidence level. 

The lack of dust rings with small radii at larger ages 
as well as the hint for an increase of the upper envelope 
of the distribution are in good accordance with the pre- 
dicted outward propagation of the dust pro duction site 
as t he result of self- or planetary stirring. Rhe e et aD 
([2007) also reported increasing radii at larger ages for a 
sample of late B- and A-type stars. The latter authors es- 
timated disk radii identically to our approach fSect. [l^ . 
thus direct comparison with our results is meaningful. 
Apart from seven very extended disks (where the radius 
estimate was based on IRAS d ata only) th e general dis- 
tribution of points in Fig. 7 of iRhee et al.l (j2007l ) is very 
si milar to our results s hown in Fig. 1161 

iKenvon fc BromlevI (|2008() predicted that the pace of 
the outward propagation in a disk depends on the disk 
mass: the more massive the disk the faster the spread 
outwards. Thus, during the active period of self-stirring 
evolution (when the expanding ring reaches the outer 
boundary of the disk), in an initally more massive disk, 
the bright ring associated with the formation of Pluto- 
sized planetesimals is located at larger radius at any 
given age. This effect offers a good explanation for the 
large scatter in dust ring radii we observe for younger 
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stars (<50Myr). In order to test this hypothesis, we es- 
timated the relative d ust masses of the disks as Mdust oc 
fdnstRl^st (see eq. 7 in lWvattl[2008l ) for the seven 30 Myr- 
old disks. We assume that the relative disk mass dis- 
tribution does not change during disk evolution, thus 
the current relative masses reflect the initial mass dis- 
tribution. For these disks we displayed the relative dust 
masses as the function of the estimated radii of the dust 
ring in Fig. [TTl The obvious trend i s fully consistent with 
predictions of iKenyon fc BromlevI (|2008| ): more massive 
d isks are located farther from the star. 

'Kc nvon Sz BromlevI (p008) derived the characteristic 
timescales of planetesimal disk formation and evolution 
processes in a disk with an initial surface density distri- 
bution of 

S - Eo(A/,)a;™(a/ao)"^/', (5) 

where So is the reference surface density at a radius of 
qq = 30 AU, while scaling factor. The refer- 

ence surface density was scaled with the stellar mass 
as So(M*) = 0.18(M,/MQ)gcm-2 (Eq = 0.18 gcm'^ 
corresponds to the minimum mass solar nebula density 
at the radius of 30 AU). According to their results, the 
timescale for the formation of the first 1000 km icy plan- 
etesimals at a radius a can be computed as: 

tiooo = 145x-/-i5(a/80AU)3(2MQ/M,)3/2[Myr]. (6) 

Adopting the radii of dust rings as the radii of the plan- 
etesimal belts and taking into account the derived ages of 
the systems, we can estimate a minimum Xm value (prac- 
tically a minimum initial surface density) that is neces- 
sary for the formation of Pluto-sized planetesimals oc- 
curring at the given radial location. Assuming a fix disk 
geometry {Rin — 10 AU, i?out = 100 AU), an estimated 
initial disk mass (Afd) and an initial disk-to-star mass ra- 
tio (Md/M^,) can also be computed (adopting the canon- 
ical gas-to-dust ratio of 100). Fig. [T51 shows the derived 
minimum Xm.min values and the minimum initial disk-to- 
star mass ratios for our sample. Eight disks (HD 3670, 
HD 15115, HD 16743, HD 30447, HD 36968, HD 113337, 
HD 170773, HD 192758) require an initial surface den- 
sity higher than that of the minimum solar mass neb- 
ula. The largest Xm.min value (3.3) was obtained for 
HD 16743. Assuming that self-stirring occurred in these 
systems and comparing the derived Md/M, ratios with 
the cumulative distribution of the corresponding ratios 
for protoplanetary disks in the Taurus and Ophiucus 
star forming regions (Fig. 7. [Andrews fc Wilha ms 2007). 
we concluded that these eight disks mi ght represent the 
high e nd of the disk mass distribution. iMustill &: Wvatti 
(|2009f) argued that a disk with a;m,min > 10 is unlikely to 
be the site of a collisonal cascade ignited by self-stirring. 
Since the highest a;ni.min hi^ our sample is 3.3, all our disks 
could be the result of self-stirring. Note, however, that 
in the calculation of dust ring radii we assumed the pres- 
ence of large blackbody grains. The presence of smaller 
grains that are ineffective emitters would have resulted 
in larger radial locations. Due to the strong dependence 
of the minimum initial surface density on the radial loca- 
tion, the uncertainty of -Rdust can affect (usually increase) 

the value of Xm.min- 

Six out of the eight disks where our calculations im- 
plied high initial surface density are younger than 60 Myr 
(see Table H]) and are located within 80 pc (Table [T]). 



These systems might be good targets for future planet 
searching surveys via direct imaging, because 1) accord- 
ing to the models, disks with high initial surface density 
are favorable for planet formation; 2) young ages make 
giant planets more easily detectable (since they are still 
bright). We note that HR8799, Fomalhaut, and ^Pic, 
systems hosting massive outer planets, harbor massive 
bright debris disks as well. 

In the present Solar System, planetesimals reside 
mainly in two spatially separated belts, the Kuiper belt 
and the main asteroid belt. Infrared and submillimeter 
observations of debris disks implied that stars harbouring 
multiple component disks rn ight be common (jChen et al.l 
120091 ISmith fc Wvattl[20lol and references therein). The 
dust replenishment processes may also be similar to those 
in our Solar System: collisions between planetesimals are 
thought to lead to the production of fresh dust in the 
outer cold rings, while besides the collisions, sublima- 
tion of icy planetesimals (scattered from the outer reser- 
voir) can also contribute to the maintenance of the inner 
warm dust belt. Several of these systems are younger 
than 100 Myr and their age coincides well with the era 
of terrestrial planet fo rmation in our Solar System (see 
lApai fc Laurettall2010l for a review), a process believed 
to be accompanied with the release of a huge amount 
of d ust grains due to collisions between large protoplan- 
ets (jKenvon fc Bromlevll2004"bl) . The warm dust around 
these young stars might eit her originate fro m terres- 
trial planet formation (e.g. ILisse "etal|[200l or from 
the collisional grinding of a rocky asteroidal belt. In 
Sect. 14.41 we found five debris systems where the SED 
can be fitted better with a two-component temperature 
model than with a single temperature one. Assuming 
that the warm component is associated with inner dust 
rings we propose that these five stars harbor two spa- 
tially separated dust rings co-located with two distinct 
planetesimal belts. All of these systems are younger than 
40 Myr. Taking into account the above mentioned argu- 
ments, these five warm dust rings might also be linked 
to the formation of terrestrial planets or collisional evo- 
lution of an asteroid belt. Since the characteristic tem- 
perature of the dust grains (150-180 K) significantly ex- 
ceeds tlie_sublimation temperature of comets (~110K, 
I Wvatti [2008D . cometary activity can also pl ay a role in 
the re plenishment of warm dust particles. I Moor et al.l 
(200§) found two young F-type members of the Tucana- 
Horologium association (HD 13246 and HD 53842) that 
harbor warm debris rings with properties (Tdust '~150- 
170 K, /dust ~0. 5-1. 7x10"'*) very similar to those of the 
currently analyzed five disks. The inner zone of the two 
Tucana-Horologium stars might also be sites of intense 
planet formation, although, no cold excess emission has 
been seen toward them. The lack of observable cold ex- 
cess may suggest a relatively dust-free outer disk perhaps 
caused by disk truncation. 

6. SUMMARY 

With the aim of investigating the properties and evo- 
lutionary trends of debris disks around F-type stars, we 
observed 82 targets with the Spitzer Space Telescope us- 
ing the IRS and MIPS instruments to obtain infrared 
spectroscopy and photometry. The core sample was se- 
lected based on hints for excess emission in earlier in- 
frared observations, and this sample was supplemented 
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by F-type members of some young kinematic groups. We 
found 27 stars that harbor debris disks, nine of which 
are new discoveries. Apart from two stars, ah of the 
observed disks exhibit excess in the IRS spectra, and 
17 disks were detected even at 160/im. In two cases, 
the emission was found to be marginahy extended on 
the 70/im MIPS images. Combining the MIPS and IRS 
measurements with additional data taken from former 
infrared space missions (IRAS, ISO) and ground-based 
submillimeter/millimeter (SCUBA, IRAM) observations, 
we compiled the spectral energy distribution of the tar- 
gets. Thanks to the large number of measurements at 
different wavelengths we achieved excellent spectral cov- 
erage for most of our debris systems. We have modeled 
the excess emission of 22 debris disks using a single tem- 
perature dust ring model and 5 debris systems with a 
two-temperature model. The latter systems may contain 
two dust rings around the star. All of our disks were 
found to be collision dominated. 

Among the 27 disks, 15 encircle stars that likely belong 
to young moving groups, offering accurate age estimates 
for these systems. We identified five new moving group 
members. HD 15745, that exhibit one of the largest 
fractional luminosities among northern debris disks, is 
a likely member of the /3Pic moving group. HD 36968, 
a system with high fractional luminosity, can be as- 
signed to th e recently discovered Octans kinematic group 
(jTorres et al. 2008), being the first member where the 
existence of a debris disk was discovered. 

In accordance with the expected trends, the fractional 
luminosity of the disks declines with time. The dis- 
tribution of data points, in particular their upper en- 
velope, seems to suggest a decay rate falling within 
the range of the model predictions (/dust, o c ^-0-3----i-0 
Dom inik & Decin 2003; Wvatt ct al. 20072 iLohne et all 
20081: iKenvon fc Bromlevlf2008t ). All disks in the sample 
seem to be cons istent with quasi ste ady-state evolution- 
ary scenarios of lWvatt etffl((2007i) . 

Assuming blackbody grains, we computed the radial 
location of the cold dust rings in all debris systems. We 
found a hint for an increase of the upper envelope of 
the radius distribution with increasing age. While older 



systems (age > 100 Myr) harbor dust rings located at 
radii of >30AU, around younger systems there are sev- 
eral disks at radial location between 10 and 30 AU. Both 
findings are in accordance wi th the predictions of s elf- or 
planetary stirring theories of IKenvon fc BromlevI ()2008f ) 
and .MustiU fc Wvatt (2009). 
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Basic properties of our target list 
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-0.13 


1 


4.25 


6480 


0.00 


N 






2 


HD 14691 


FOV 


5.43 


29.7 


-0.12 


1,2 


4.25 


6800 


0.00 


N 






1 


HD 15060 


F5V 


7.02 


76.0 


-0.14 


1 


4.00 


6260 


0.00 


N 






1 


HD 15115 


F2 


6.79 


45.2 


-0.06 


1 


4.25 


6780 


0.00 


N 






1 


HD 15745 


FO 


7.47 


63.5 


-0.10 


1 


4.25 


6860 


0.00 


N 






1 


HD 16743 


F0/F2III/IV 


6.78 


58.9 


-0.13 


1 


4.25 


7000 


0.00 


Y 


216.2 




1 


HD 17390 


F3IV/V 


6.48 


48.0 


0.03 


1 


4.25 


6840 


0.00 


N 






1 


BD+49° 896 


F4V 


9.68 


(175.0) 


-0.15 


3 


4.25 


6740 


0.06 


N 






1 


HD 20759 


F5V 


7.70 


76.8 


-0.48 


1 


4.00 


6280 


0.00 


N 






1 


HD 24636 


F3IV/V 


7.13 


54.1 


-0.11 


1 


4.25 


6820 


0.00 


N 






2 


HD 25570 


F2V 


5.45 


34.9 


-0.29 


1 


4.00 


6760 


0.00 


N 






1 


HD 25953 


F5 


7.83 


55.2 


-0.22 


1 


4.25 


6240 


0.00 


N 






2 


HD 27429 


F3:V... 


6.11 


48.3 


0.01 


1 


4.00 


6720 


0.00 


Y 






1 


HD 30447 


F3V 


7.85 


80.3 


-0.21 


1 


4.25 


6800 


0.00 


N 






1,2 


HD 32195 


F7V 


8.14 


61.0 


-0.13 


1 


4.25 


6180 


0.00 


N 






2 


HD 30743 


F3/F5V 


6.27 


33.8 


-0.38 


1,2,4,5 


4.25 


6440 


0.00 


N 






1 


HD 33081 


F7V 


7.04 


50.6 


-0.19 


1 


4.25 


6360 


0.00 


N 






1 


HD 33276 


F2IV 


4.81 


165.0 


0.27 


6,7 


3.25 


6920 


0.00 


Y 


0.300 




1 


HD 34739 


F7IV/V 


9.33 


(121.0) 




4.25 


6280 


0.00 


N 






1 


HD 35114 


F6V 


7.39 


48.3 


-0.18 


1 


4.25 


6200 


0.00 


N 






2 


HD 35841 


F5V 


8.91 


( 96.0) 






4.25 


6460 


0.00 


N 






1,2 


HD 36248 


F8 


8.05 


73.6 


0.11 


1 


4.25 


5960 


0.06 


Y 


3.800 




1 


HD 36968 


F2V 


9.02 


(140.0) 






4.25 


6880 


0.00 


N 






1,2 


HD 37402 


F6V 


8.38 


78.2 


-0.15 


1 


4.25 


6160 


0.00 


N 






2 


HD 38905 


F6/F7V 


9.73 


(140.0) 






4.25 


6240 


0.04 


N 






1 


HD 47412 


F2 


6.82 


111.0 


-0.03 


1 


4.00 


6300 


0.00 


Y 


0.500 




1 


HD 48391 


F5 


7.89 


58.3 


-0.20 


1 


4.25 


6160 


0.00 


N 






2 


HD 50571 


F7IH-IV 


6.11 


33.6 


-0.02 


1,2 


4.25 


6480 


0.00 


N 






1 


HD 55003 


F2 


7.04 


60.1 


-0.02 


1 


4.00 


6340 


0.00 


Y 




2 


1 


HD 56099 


F8 


7.62 


86.7 


-0.03 


1 


4.00 


6060 


0.00 


Y 


0.130 


4 


1 


HD 58853 


F5V 


9.07 


122.4 






4.25 


6280 


0.00 


N 






1 


HD 61518 


F5V 


7.88 


62.0 


-0.19 


1 


4.25 


6340 


0.00 


N 






2 


HD 67587 


F8 


6.65 


46.8 


-0.19 


1 


4.00 


5980 


0.00 


Y 




2 


1 


HD 69351 


F8 


7.16 


80.0 


0.08 


1 


4.00 


5980 


0.00 


Y 


1.400 


3 


1 


HD 79873 


F5 


6.73 


68.3 


0.13 


1 


4.00 


6400 


0.00 


Y 


2.100 


3 


1 


HD 82821 


F8 


8.69 


74.7 






4.25 


6140 


0.00 


Y 




2 


1 


HD 86146 


F6Vs 


5.11 


28.1 


0.02 


1 


4.00 


6360 


0.00 


Y 




5 


1 


PPM 7774 


F5 


8.96 


(116.0) 






4.25 


6520 


0.00 


N 






1 


HD 103257 


F2V 


6.62 


63.4 


-0.39 


1 


4.00 


6980 


0.00 


N 






1 


HD 107067 


F8... 


8.69 


66.0 


-0.09 


1,8,9 


4.25 


6060 


0.00 


N 






1 


HD 108102 


F8... 


8.12 


95.1 


-0.13 


1 


4.00 


6060 


0.00 


N 






1 


HD 113337 


F6V 


6.01 


36.9 


0.06 


1,10 


4.25 


6600 


0.00 


Y 


119.7 




1 


HD 114905 


F7V 


6.83 


61.4 


-0.21 


1 


4.00 


6280 


0.00 


N 






1 


HD 117360 


F6V 


6.52 


35.2 


-0.34 


1,2 


4.25 


6400 


0.00 


Y 


22.40 




1 


HD 120160 


FOIV/V 


7.67 


138.3 


0.02 


1 


3.75 


6760 


0.00 


N 






1 


HD 122106 


F8V 


6.36 


77.5 


0.08 


1,11 


3.75 


6280 


0.00 


N 






1 


HD 122510 


F6V 


6.18 


38.2 


-0.11 


1,2 


4.25 


6600 


0.00 


Y 


1.900 




1 


HD 124988 


FO 


6.88 


96.0 


-0.09 


1 


4.00 


6880 


0.00 


N 






1 


HD 125451 


F5IV 


5.41 


26.1 


-0.00 


1,12 


4.25 


6660 


0.00 


N 






1 


HD 127821 


F4IV 


6.10 


31.8 


-0.18 


1 


4.25 


6660 


0.00 


N 






1 


HD 131495 


F2 


6.87 


72.3 


-0.13 


1 


4.00 


6460 


0.00 


N 






1 


HD 134150 


F8 


9.84 


(166.0) 






4.25 


6460 


0.03 


N 






1 


HD 136580 


F5 


6.90 


41.0 


-0.15 


1,13 


4.25 


6180 


0.00 


Y 




2 


1 


HD 136407 


F2V 


6.14 


56.7 


-0.08 


1 


4.00 


6640 


0.00 


Y 


44.40 


3 


1 


HD 138100 


FO 


6.69 


57.8 


-0.07 


1 


4.00 


6720 


0.00 


Y 




2 


1 


HD 139798 


F2V 


5.76 


35.7 


-0.22 


1,12 


4.00 


6700 


0.00 


Y 




2 


1 


HD 143840 


FIV 


8.11 


132.3 


-0.03 


1 


4.00 


6680 


0.47 


N 






1 


HD 145371 


GO 


9.46 


(144.0) 






4.25 


6480 


0.00 


N 






1 


HD 151044 


F8V 


6.48 


29.3 


-0.02 


1,9,12,13,14,15 


4.25 


6060 


0.00 


N 






1 


HD 153377 


F2 


7.55 


64.0 


-0.16 


1 


4.25 


6620 


0.00 


N 






1 


HD 155990 


F8 


8.07 


57.7 


-0.18 


1 


4.25 


6100 


0.00 


N 






2 


HD 156635 


F8 


6.66 


40.3 


-0.09 


1,16 


4.25 


6160 


0.00 


Y 






1 


HD 170773 


F5V 


6.22 


37.0 


-0.05 


1,2 


4.25 


6640 


0.00 


N 






1 


HD 184169 


F2 


8.20 


(83.0) 






4.25 


6540 


0.00 


N 






1 


HD 183577 


F6V 


6.48 


41.6 


-0.24 


1 


4.00 


6100 


0.00 


Y 






1 


HD 185053 


F5/F6V 


8.83 


(67.0) 






4.25 


5960 


0.24 


N 






1 


HD 189207 


F2 


8.08 


116.3 


-0.06 


1 


4.00 


6780 


0.00 


N 






1 


HD 192486 


F2V 


6.55 


44.7 


-0.18 


1 


4.25 


6820 


0.00 


N 






1 


HD 192758 


FOV 


7.03 


(62.0) 


-0.06 


17 


4.25 


7080 


0.00 


N 






1 


HD 195952 


F3V 


8.12 


153.4 


-0.06 


1 


3.75 


6240 


0.06 


N 






1 


HD 199391 


F0/F2IV 


7.12 


78.7 






4.00 


7160 


0.00 


Y 


4.800 




1 


PPM 171537 


F8 


9.23 


(90.0) 






4.25 


5960 


0.00 


N 






1 


HD 204942 


F7V 


8.23 


84.9 


-0.14 


1 


4.25 


6220 


0.00 


N 






1 
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TABLE 1 — Continued 







(2) 


(3) 


(4) 


(5) 




C7) 


V"/ 


("9) 


dol 




C12~I 


(13) 




ID 


SpT 


V 


D 


fFe/Hl 


Ref. 


lop" P" 


T a 
^ erf 


-n V 


Mult. 


Sep. 


Ref. 


Sel. cr. 








[llldgj 
















I 1 






HD 


205674 


F3/F5IV 


7.19 


51.8 


-0.23 


1 


4.25 


6780 


0.00 


N 






1 


HD 


206554 


F5 


7.12 


65.9 


-0.23 


1 


4.00 


6400 


0.00 


N 






1 


HD 


206893 


F5V 


6.69 


38.3 


-0.06 


1,18 


4.25 


6520 


0.00 


N 






1 


HD 


207889 


F5 


7.20 


49.6 


-0.11 


1 


4.25 


6520 


0.00 


N 






1 


HD 


210210 


FlIV 


6.08 


88.2 






3.75 


7100 


0.00 


N 






1 


HD 


213429 


F7V 


6.15 


25.4 


-0.01 


1,19 


4.25 


6040 


0.00 


Y 


0.072 




1 


HD 


213617 


FIV 


6.43 


50.3 


-0.11 


1 


4.25 


7020 


0.00 


N 






1 


HD 


218980 


F 


8.58 


(105.0) 






4.25 


7060 


0.43 


N 






1 


HD 


221853 


FO 


7.35 


68.4 


-0.05 


1,18 


4.25 


6760 


0.00 


N 






1 



Note. — Col.(l): Identification. Col. (2): Spectr al type. Col. (3): V magnitude. Col. (4): Distance. Parentheses indicate photometric 
distances, otherwise Hipparcos distances from I vairL ccuwcn (2007) are used. Col. (5): Metallicity. Literature data are used; if more t han 
one observation is available the average of the [Fe/H] is quoted. Col. (6): References for metallicity data. DI Holmberg et al, | l|200^, 2 ) 
[Grav ct al. (2006), 3) Bocsgaard et al. (1988b), 4) Takcda (2007), 5) Edvardsson ct al. (1993), 6) Bcrthct (1990). 7) Ccnarro ct a.1. (20071), 
8) Bocsgaard (1987), 9) Fricl & Bocsgaard (1992), 10) Bocsgaard & Tripicco (1986), 11) Balachandran (1990), 12) Bocsgaard ct at] 
(1988311 . 131 iValenti fc F ischer (2005), 14) Monicr (2005), 15) IBoesgaard fc Frieli 1 119901') . 16) iReddv et al.l J2003i), 17) Me tallicity was 
derived based on uvby ~^ a,uck & Mcrmilliod 1997) photometric data, using the calibration described by [Holmberg et al.l (2007). 18) 
ISatt'e et al.l (120081 ). 19) IBoesgaard et al.i (|2004 ). Col. (7): Surface gravity values fixed in the course of fitting stellar atmospheric models. 
Col. (8): Derived effective temperature. Col. (9): Interstellar extinction. Col. (10): Multiplicity. Col. (ll): Separa tion of t he compo- 
nents if the object is in multiple system. Col. (12): References for multiplicity data. 1) See Sect. |431 2) Frankow ski et aLl ((20071), 3) 
Ipommangct & Nvs (2002), 4) Balcga et al. (2007), 5) Batten & Morbcv (1980), 6) Pourbaix (2000). Col. (13): Selection criteria: 1) 
the star was suspected to display IR excess based on earlier IR observations; 2) the star was selected because of its kinematic group 
membership. 
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TABLE 2 

MIPS PHOTOMETRY 



S/^iiTT'fi TTi 


AOR KEY 




24/j.rQ 






70^m 






1 6n/(m 

A-KjyjLLlll. 




J. I 


r 24 L^i-'-'Jyj 




X24 


r 70 [iiijyj 






17 1 f/^ \yYt Tifl 

160 [-liA'jyj 




X160 


HD 3670 


15010816 


12.8±0.5 


10.4 


3.9 


134.9±10.4 


1.1 


12.8 


77.2±13.4 


0.2 


5.7 




HD 14691 


14996736 


114.4±4.5 


116.9 


-0.4 


17.7±5.0 


12.7 


0.9 










HD 15060 


23050496 


34.6±1.3 


34.6 


-0.0 


14.8±3.5 


3.7 


3.1 










HD 15115 


10885888 


58.3±2.3 


33.5 


9.7 


451.9±32.6 


3.6 


13.7 


217.3±27.8 


0.7 


7.7 




HD 15745 


10886400 


169.4±6.7 


17.3 


22.3 


741.0±52.6 


1.8 


14.0 


230.8±29.9 


0.3 


7.6 




HD 16743 


15002624 


50.3±2.0 


30.1 


9.1 


368.8±26.5 


3.2 


13.7 


174.7±24.5 


0.6 


7.0 


5 


HD 17390 


10887168 


45.3±1.8 


41.8 


1.5 


255.2±18.9 


4.5 


13.2 


210.0±28.3 


0.9 


7.3 




BD+49° 896 


15008000 


2.4±0.1 


2.6 


-0.5 


0.5±4.9 


0.2 


0.0 










HD 20759 


23051008 


19.0±0.7 


18.6 


0.4 


9.6±4.0 


2.0 


1.8 










HD 24636 


23051520 


42.3±1.6 


24.3 


9.7 


35.1±4.7 


2.6 


6.8 










HD 25570 


15002880 


115.6±4.6 


113.7 


0.3 


170.5±13.3 


12.4 


11.8 


115.0±28.4 


2.5 


3.9 


6 


HD 25953 


15009792 


16.7±0.6 


16.4 


0.3 


-5.0±4.4 


1.8 


-1.5 


-6.3±41.3 


0.3 


-0.1 




HD 27429 


14996992 


63.8±2.5 


61.7 


0.6 


6.8±6.9 


6.7 


0.0 










HD 30447 


10887680 


30.1±1.2 


12.5 


13.8 


289.8±21.1 


1.3 


13.6 


120.3±17.6 


0.2 


6.8 




HD 32195 


15009536 


16.4±0.6 


12.6 


4.9 


17.0±3.9 


1.3 


3.9 


-21.3±20.3 


0.2 


-1.0 




HD 30743 


15003136 


66.6±2.6 


65.7 


0.2 


13.9±7.5 


7.2 


0.8 


5.4±17.2 


1.4 


0.2 


5,6 


HD 33081 


14997248 


33.8±1.3 


32.6 


0.6 


35.8±5.1 


3.5 


6.2 








HD 33276 


15003392 


198.9±7.9 


181.2 


1.8 


26.8±9.2 


19.6 


0.7 


50.9±42.5 


4.0 


1.1 


6 


HD 34739 


10888192 


4.0±0.1 


4.1 


-0.4 


-1.7±4.9 


0.4 


-0.4 








2,3,4 


HD 35114 


11260160 


32.2±1.2 


25.2 


4.6 


20.2±4.2 


2.7 


4.0 










HD 35841 


15011840 


18.4±0.7 


5.4 


17.0 


172.1±13.6 


0.5 


12.5 


37.5±14.2 


0.1 


2.6 




HD 36248 


14997504 


16.8±0.6 


16.3 


0.4 














1 


HD 36968 


15012096 


8.7±0.3 


4.1 


12.1 


148.1±12.5 


0.4 


11.7 


105.0±16.6 


0.09 


6.3 


5 


HD 37402 


15011072 


10.4±0.4 


10.3 


0.1 


-1.6±4.2 


1.1 


-0.6 


-19.5±12.0 


0.2 


-1.6 


5 


HD 38905 


15012352 


3.1±0.1 


2.9 


1.3 


-0.2±7.4 


0.3 


-0.0 








2,3,4 


HD 47412 


14997760 


38.3±1.5 


39.5 


-0.5 


5.8±5.4 


4.3 


0.2 










HD 48391 


15010048 


15.6±0.6 


16.0 


-0.4 


17.1±5.6 


1.7 


2.7 


-26.9±20.5 


0.3 


-1.3 




HD 50571 


15003904 


70.4±2.8 


69.4 


0.2 


248.8±18.7 


7.5 


12.8 


214.6±36.0 


1.5 


5.9 


6,7 


HD 55003 


23054848 


31.6±1.2 


32.9 


-0.7 


2.5±3.6 


3.5 


-0.2 








HD 56099 


15004160 


22.3±0.8 


21.8 


0.4 


10.6±7.6 


2.3 


1.0 








3,4 


HD 58853 


15004416 


5.0±0.2 


5.1 


-0.2 


6.6±6.5 


0.5 


0.9 








3,4 


HD 61518 


15010304 


15.4±0.6 


15.1 


0.4 


-9.0±5.4 


1.6 


-1.9 


3.9±17.2 


0.3 


0.2 


HD 67587 


14998016 


56.3±2.2 


58.3 


-0.6 


12.9±9.6 


6.3 


0.6 










HD 69351 


14998272 


31.7±1.2 


32.8 


-0.6 


-1.6±5.9 


3.5 


-0.8 










HD 79873 


14998528 


40.7±1.6 


41.2 


-0.2 


-2.1±4.8 


4.4 


-1.3 










HD 82821 


15004672 


7.5±0.3 


7.8 


-0.5 


4.0±8.0 


0.8 


0.3 


-4.5±13.3 


0.1 


-0.3 


5 


HD 86146 


14998784 


183.3±7.3 


186.2 


-0.3 


21.4±4.6 


20.2 


0.2 










PPM 7774 


15011584 


4.9±0.2 


5.1 


-0.6 


-0.3±7.3 


0.5 


-0.1 


-90.6±43.5 


0.1 


-2.0 


5 


HD 103257 


14999040 


36.7±1.4 


35.3 


0.7 


12.7±7.3 


3.8 


1.2 










HD 107067 


15008768 


8.2±0.3 


8.4 


-0.6 


-5.4±10.6 


0.9 


-0.6 


-21.5±10.4 


0.1 


-2.0 




HD 108102 


15009024 


14.2±0.5 


13.9 


0.3 


-0.0±7.0 


1.5 


-0.2 


21.9±11.6 


0.3 


1.8 




HD 113337 


14999296 


74.7±2.9 


71.3 


0.9 


178.2±13.3 


7.7 


12.7 










HD 114905 


14999552 


40.3±1.6 


40.9 


-0.3 


3.1±4.1 


4.4 


-0.3 










HD 117360 


19890432 


51.0±2.0 


55.4 


-1.6 


6.6±6.0 


6.0 


0.0 










HD 120160 


15004928 


14.6±0.6 


14.1 


0.7 


30.7±7.3 


1.5 


3.9 


1.2±22.7 


0.3 


0.0 




HD 122106 


15005184 


58.8±2.3 


59.2 


-0.1 


1.8±7.7 


6.4 


-0.5 








4 


HD 122510 


15000064 


70.6±2.8 


70.7 


-0.0 


8.8±5.3 


7.7 


0.2 










HD 124988 


23062016 


28.8±1.1 


28.3 


0.3 


9.7±4.1 


3.0 


1.5 










HD 125451 


15008256 


126.2±5.0 


121.9 


0.6 


64.6±8.1 


13.2 


6.3 


-2.0±11.2 


2.7 


-0.4 


6 


HD 127821 


15005440 


69.8±2.7 


69.9 


-0.0 


360.3±25.9 


7.6 


13.5 


295.5±37.2 


1.5 


7.8 


6 


HD 131495 


15000320 


34.1±1.3 


35.4 


-0.7 


-4.8±4.8 


3.8 


-1.8 










HD 134150 


10888960 


2.4±0.1 


2.4 


0.3 


-6.5±7.4 


0.2 


-0.9 


3.5±15.7 


6.05 


0.2 


5 


HD 136580 


15000576 


40.1±1.6 


40.4 


-0.1 


-0.9±4.2 


4.4 


-1.2 










HD 136407 


15005696 


66.8±2.6 


64.0 


0.8 


10.3±7.9 


6.9 


0.4 


-6.0±22.8 


1.4 


-0.3 


6 


HD 138100 


15000832 


36.5±1.4 


37.5 


-0.5 


0.8±4.2 


4.0 


-0.7 










HD 139798 


15008512 


88.9±3.5 


89.9 


-0.2 


10.7±5.9 


9.8 


0.1 


-20.7±10.6 


2.0 


-2.1 


6 


HD 143840 


10889472 


15.4±0.6 


14.7 


0.9 














1 


HD 145371 


15007744 


3.3±0.1 


3.2 


0.4 


0.1±7.4 


0.3 


-0.0 








2,3,4 


HD 151044 


10890240 


68.3±2.7 


63.7 


1.3 


95.6±7.9 


6.9 


11.2 


47.4±11.7 


1.4 


3.9 


6 


HD 153377 


15005952 


18.5±0.7 


18.8 


-0.3 


-2.0±5.6 


2.0 


-0.7 


26.9±24.0 


0.4 


1.1 




HD 155990 


15010560 


14.8±0.6 


14.5 


0.4 


5.7±5.1 


1.5 


0.8 




0.3 


2.3 


4 


HD 156635 


15006208 


48.2±1.9 


50.1 


-0.7 


6.9±6.1 


5.4 


0.2 


56.7±75.4 


1.1 


0.7 




HD 170773 


10890752 


65.3±2.6 


60.5 


1.5 


787.9±56.0 


6.6 


13.9 


692.3±83.8 


1.3 


8.2 


6,7 


HD 184169 


10891264 


10.5±0.4 


10.3 


0.2 


-12.1±12.4 


1.1 


-1.0 








2,3,4 


HD 183577 


15001600 


63.0±2.5 


61.7 


0.4 


4.8±4.6 


6.7 


-0.4 










HD 185053 


10892032 


22.8±0.9 


9.8 


13.4 














1 


HD 189207 


15006464 


10.1±0.4 


9.8 


0.5 


5.9±5.4 


1.0 


0.8 


-46.4±30.6 


0.2 


-1.5 




HD 192486 


15001856 


43.3±1.7 


42.1 


0.5 


-5.2±4.6 


4.6 


-2.1 










HD 192758 


10892800 


42.1±1.6 


23.4 


10.1 


452.4±32.5 


2.5 


13.8 


200.7±26.3 


0.5 


7.5 




HD 195952 


15006720 


12.9±0.5 


13.0 


-0.2 


1.4±8.1 


1.4 


0.0 








3,4 


HD 199391 


15006976 


22.3±0.9 


24.4 


-1.7 


-7.5±10.2 


2.6 


-0.9 


130.7±27.8 


0.5 


4.6 


5 


PPM 171537 


15011328 


5.5±0.2 


5.3 


0.6 


-2.4±7.7 


0.5 


-0.3 








3,4 


HD 204942 


15007232 


10.9±0.7 


11.5 


-0.6 


-13.8±11.0 


1.2 


-1.3 








3,4 


HD 205674 


15007488 


31.2±1.2 


24.1 


4.8 


232.5±17.9 


2.6 


12.8 


185.6±26.4 


0.5 


6.9 
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TABLE 2 — Continued 



^rtiirrT" TTi 


A OR KEY 




24/j.ni 












1 fiO / ( tn 

_1_ U vy iXlll 






r 24 L^i-'-'Jyj 




X24 


r 70 [iiijyj 




X70 


^ 160 i^iijyj 




X160 


HD 206554 


23067648 


29.0±1.1 


29.2 


-0.1 


1.4±6.3 


3.1 


-0.2 








3 


HD 206893 


10893312 


44.2±1.7 


40.6 


1.6 


265.7±20.0 


4.4 


13.0 


193.1±26.2 


0.9 


7.3 




HD 207889 


23068160 


27.0±1.0 


25.5 


1.1 


4.7±3.8 


2.7 


0.4 








3 


HD 210210 


15002112 


53.2±2.1 


52.9 


0.1 


5.0±4.0 


5.7 


-0.1 










HD 213429 


21840896 


92.2±3.6 


86.3 


1.3 


22.2±4.1 


9.4 


3.0 










HD 213617 


10893824 


44.3±1.7 


41.8 


1.1 


119.2±10.0 


4.5 


11.3 


96.4±22.9 


0.9 


4.1 


5 


HD 218980 


13238784 


11.3±0.4 


8.2 


5.9 














1 


HD 221853 


10894848 


78.5±3.2 


20.3 


17.7 


336.5±24.7 


2.2 


13.5 


105.3±20.4 


0.4 


5.1 





Note. — Col.(l): Identification. Col. (2): AOR Key for MIPS measurement. Using the AOR key one can query additional details for each observation (e.g. measurement 
setups) from the Spitzer Data Archive at the Spitzer Science Center. Col. (3-11). Measured and predicted flux densities with their uncertainties and the significance of 
the excesses at 24/70/160/^m. The quoted uncertainties include the calibration uncertainties as well. Col. (12): Notes. 1: Nebulosity. At 24/^m a small aperture with 3.5" 
radius was used. At 70/^m and 160/J.m no photometric values were quoted; 2; A bright nearby source contaminated the aperture photometry at 24^im — PSF photometry 
was used; 3; Bright nearby sources contaminated the photometry of our target at 70ftm. In order to remove the contribution of these background objects, we fitted PSF 
to these sources and subtracted their emission before the photometry for our target was performed. 4: Nearby bright object at 160/^m in the aperture - no photometry 
is given; 5: A bright nearby source contaminated the annulus at 160/im. The background source was masked out in the course of photometry; 6: The ghost image was 
subtracted at 160/im before the photometry was performed. 7: The source is marginally extended at 70/im and the photometry was derived using the fitted profile (see 
Sect. f43t . 
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TABLE 3 
Additional photometric data 



Source ID 


Instrument 


Wavelength 








[fim] 


[mjy] 


HD 151044 


ISO/ISOPHOT 


60 


111.0±13.0 


HD 151044 


ISO/ISOPHOT 


90 


lOl.OilO.O 


HD 170773 


ISO/ISOPHOT 


60 


570.0±35.0 


HD 170773 


ISO/ISOPHOT 


90 


771.0±54.0 


HD 213617 


ISO/ISOPHOT 


60 


110.0±13.0 


HD 213617 


ISO/ISOPHOT 


90 


121.0±14.0 


HD 15745 


IRAM/MAMB02 


1200 


1.3±0.6 


HD 25570 


IRAM/MAMB02 


1200 


0.6±0.6 


HD 113337 


IRAM/MAMB02 


1200 


0.4±0.3 



TABLE 4 

Additional spectroscopic information and derived kinematic 
properties 



Source ID 


Vr 


U 


V 


W 


vsini 


EWl, 




[kms-1] 


[kms-i] 


[kms-i] 


[kms^-*-] 


[kms^-^] 


[mA] 


HD 15745 


2.5±3.3 


-10.4±2.5 


-15.3±1.9 


-7.9±1.2 


50 




HD 16699 


16.2±0.2 


-23.5±1.5 


-15.0±0.4 


-10.3±0.3 


20 


40±10 


SAO 232842 


15.8±0.3 










220±10 


HD 16743 


17.3±8.5 


-23.3±0.5 


-15.4±4.6 


-11.3±7.2 


100 




HD 24636 


14.0±0.7 


-8.8±0.3 


-20.5±0.6 


-1.6±0.5 


30 


50±10 


HD 36968=^ 


15.0±2.0 


-14.8±2.1 


-6.6±2.0 


-8.5±1.3 






HD 170773 


-17.5±1.5 


-22.2±1.5 


-4.6±0.2 


-15.0±0.5 


50 




HD 205674 


-1.4±1.3 


-3.0±0.8 


-25.2±1.1 


-14.6±1.1 


30 


20±10 


HD 206893 


-11.8±1.6 


-19.2±0.9 


-7.2±0.8 


-2.7±1.1 


33 




HD 221853 


-8.4±1.8 


-11.9±0.6 


-22.4±1.4 


-7.3±1.5 


40 





Note. — Col.(l): Identification. Col. (2): Heliocentric radial velocity. Col. (3-5): Galactic space velocity components of the star. In the 
calculation of the Galactic space velocity we used a right-handed coordinate system {U is positive towards the Galactic centre, V is positive in the 
direction of galactic rotation and W is positive towards the galactic North pole) and followed the general recipe described in "The Hipparcos and 
Tycho Catalogues" (ESA 1997). Col. (6) Projected rotational velocity of the star. Col. (7) Measured lithium equivalent width. 

^ HD 36968 was measured in the framework of a previous program that carried out with the 2.3 m ANU-telescope at the Siding Spring Observatory 
(Australia), using the Double Beam Spectrograph (for details see IMoor et al.|[2006l) . 
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TABLE r> 

Synthetic IRS photometry for stars with infrared excess 



Source ID 


AOR KEY 


8-lOMm 


10- 


-12/j,m 


12- 


-14^im 


14-16Mm 


Sy nt hct i 
16-18/^m 


c ptiotomctry" 
18-20/^m 


[mJy] 
20-23/^ m 


23-26 /^m 


26- 


-29^i] 




29- 






32-35^111 


HD 


3670 


15025664 


71.8±9.0 


48 


.9 ±5,0 


35. 8 ±4. 4 


26.7±2.6 


20,1±2,0 


16,6±2,5 


14,7±2,0 


13.2±1.0 


14, 


.1±1, 


.3 


15, 


.9±1.4 


20 


,9±2.6 


HD 


15060 


23050752 


227.6±30.4 


149 


.9±15.9 


109.9±10.0 


89.4±7.5 


69,5±6,1 


53,2±6,6 


43,9±5,8 


33.3±2.4 


26. 


.5±2, 


.1 


21. 


.3±2. 




17 


,5±3.7 


HD 


15115 


10SS5632 


218.0±36.8 


148 


.7±15.6 


111, 


.2±10.0 


87.0±7.7 


70,8±4,6 


61,0±2,3 


53,7±2,7 


54.9±2.5 


60. 


.3±2, 


.8 


71. 


.8±5. 


1 


93 


,5±9.8 


HD 


15745 


10SS6144 


123.5±15.8 


83 


.9 ±8, 3 


65.2±4.1 


67.7±3.9 


81,2±7,5 


110,2±9,7 


143,2±19,2 


218,9±16,7 


276. 


.8±19.4 


335. 


.2±22,7 


431 


,0±30.5 


HD 


16743 


15018240 


208.5±29.0 


142 


.5±15,0 


105 


,9±7,3 


79,1±6,0 


62,9±3,8 


54,5±2,9 


50,7±2,0 


50,4±1,8 


52, 






57. 






70 


,4±5.5 


HD 


17390 


10886912 


278.2±30.6 


196 


.6±21,9 


141. 


6 ±12,9 


103,7±8,2 


80,5±5,3 


63,3±4,5 


52,2±4,2 


42, 8 ±2, 6 


37. 


^±1 




36. 


2±"l 


9 


38 


,7±3,0 


HD 


24636 


23051776 


163.1±19.9 


111 


.5±11,1 


84,2±7,8 


68,7±4,3 


57,4±3,0 


50,2±3,3 


46,6±3,1 


44,2±1,4 


43. 


.6±1. 


2 


44. 


2±1. 


.9 


44 


,0±4,3 


HD 


25570 


15018496 


755.7±143.7 


530 


.2 ±56. 6 


383. 


5 ±36, 7 


285,3±21,6 


222.8±16.0 


178,1±12,9 


142,9±13,5 


111,2±7,8 


92. 


.1±4. 


6 


79. 


0±6. 


3 


72 


,2 ±6, 6 


HD 


30447 


10887424 


83.0±11.1 


55 


.9±5.8 


41, 


.7±4,3 


34,0 ±2, 5 


29.1±2,0 


26,4±2,5 


28,5±3,0 


35,7±1,9 


45. 


.2±3, 


.9 


58. 


5 ±5. 


4 


80 


,4±7,2 


HD 


32195 


15024384 


86.6±10.7 


58 


.8±6.7 


43,4±4,2 


3 1,5 ±2, 3 


25.3±2,3 


20,6±2,3 


17,4±2,1 


15, 5 ±1.1 


13, 


.5±1. 


.1 


11. 


.9±1. 


1 


12 


,0±2,0 


HD 


33081 


15013120 


219.5±28.6 


151 


.5±16.1 


109,3±9,8 


82,0±5.9 


64,3±4,3 


51,6±5,4 


41,7±3,8 


33.5±2.1 


27, 


.6±1, 


.8 


23, 


.5±1. 


.8 


20 


,2±1,S 


HD 


35114 


26362112 


171.5±25.5 


110 


.1±13.8 


81, 


.5±8.4 


62. 3 ±5, 9 


50,7±3,3 


40.8±4.1 


35.9±3,5 


29.9±1.7 


26. 


.5±1, 


.9 


24. 


.1±2. 


.1 


22 


,5±4,7 


HD 


35841 


15026688 


36.8±5.0 


24 


.7±3.0 


17, 


.7±1.8 


13,7±1,3 


11,6±1,3 


10.7±1.7 


13.4±2,3 


20.7±3.0 


28, 


.4±3, 


.1 


37, 


.1±3. 


.3 


50 


,7±5,4 


HD 


36968 


15026944 


28.2±3.4 


19 


.4±2.2 


14. 


.5±1.7 


9,9±1,4 


9,4±1,3 


7,3±1,4 


7,0±1,7 


9,3±1,0 


12, 


.3±1. 


.6 


17, 


.2±2, 


.4 


24 


,9 ±3, 9 


HD 


50571 


15019520 


422.2±28.8 


311 


.6±34.2 


225. 


.3±21,0 


172.7±12.3 


134.0±10.0 


104,5±7,0 


84,5±6,0 


68, 5 ±5,0 


56. 


.6±2, 




49. 


.1±3. 


.3 


50 


,7±4.9 


HD 


113337 


15015168 


475.6±64.0 


321 


.6 ±35. 7 


232. 


.5 ±2 1.1 


176.5±12.8 


136,9±9.6 


108,5±6,6 


87,5±7,7 


72, 8 ±4, 5 


62, 


.4±1. 


.9 


57. 


.4±1. 


.6 


58 


,3±3,8 


HD 


120160 


15020544 


94.9±13.4 


64 


.1±6,8 


46. 


0±5.2 


34.8±3.0 


27.4±2,1 


22,6±2,1 


17,9±2,2 


14, 9 ±1,3 


12. 


.7±0. 


.9 


11. 


.5±0, 


.9 


12 


.1±1.4 


HD 


125451 


15027968 












297.1±18.8 


235.7±15.8 


187.1±12,5 


152,4±12,0 


121,0±6,3 


98, 


.9±5. 


.8 


82. 


.0±5, 


.4 


73 


,2±4.0 


HD 


127821 


15020800 


475. 5± 63. 2 


319 


.8 ±34, 7 


230. 


8±21.9 


176.1±13.0 


136,6±9.7 


108,9±7,6 


88,6±7,5 


7 1,3 ±3, 5 


59, 


.3±2, 




53, 


.9±2. 


.6 


52 


,6±4.3 


HD 


151044 


10889984 


411.6±58.7 


277 


.6±26.0 


209. 


1±15.4 


156.4±15.3 


120,2±9,3 


95,5±4,5 


77,5±4,7 


64,7±4.1 


54, 


.6±2, 




47. 


.8±2. 


.4 


47 


,0±2.7 


HD 


170773 


10890496 


406. 5± 53. 6 


276 


.6 ±30.0 


200. 


.1±19.1 


149.5±12.4 


117,4±9,4 


94,5±6,3 


77,5±6,4 


63.6±3.2 


55, 


.9±1. 


.8 


54. 


.6±1. 


.8 


62 


,1±3,4 


HD 


192758 


10892544 


157.0±20.4 


108 


. 5 ± 1 1 . 2 


79, 


.6±7.0 


61.7±5.1 


49,5±3,3 


42,7±2,3 


39,8±1,9 


43.0±2.4 


49, 


.1±3, 


.1 


61. 


4±4. 


8 


81 


,9±7,2 


HD 


205674 


15022592 


162. 3± 2 1.3 


110 


.1±11.7 


80. 


.6±7.1 


61.0±6.2 


47,5±4,0 


37,4±3,6 


32,9±2,8 


31.0±1.3 


30. 


.4±1. 


.8 


34. 


5 ±2. 


6 


39 


,9±2,7 


HD 


206893 


10893056 


277.5±36.2 


188 


.5±20.8 


135, 


.6±13.0 


102.9±8.6 


80,2±5,9 


63,3±4,9 


51,7±5,3 


42.7±2.0 


38. 


.2±1, 


.3 


38. 


0±0. 


9 


42 


,2±2,4 


HD 


213429 


15017984 


575.4±84.9 


391 


.4±42.8 


284, 


.5±23.8 


218.8±18.7 


166.9±13.2 


130,1±11,4 


105.0±9,3 


78.4±5.7 


61. 


.6±5, 




48. 


.7±4. 




38 


,3±4,9 


HD 


213617 


10893568 


285.9±38.2 


191 


.0±21.0 


137. 


.9±12.8 


102.9±8.3 


80.9±5.0 


66.6±5.4 


53.2±4.2 


43.8±2.2 


38. 


.2±1. 


3 


36. 


.5±1. 


3 


36 


.3±3.4 


HD 


221853 


10894592 


128.5±18.0 


87. 


.1±9.9 


63. 


,3±4.9 


53.5±3.3 


48.4±1.3 


50.9±3.0 


64.3±8.8 


92.9±9.8 


119 


|.6±8.7 


140 


i.l±6.3 


170 


.9±12.2 



No 1 l-J.. — Col , ( 1 ) : Identification , Col , (2) : AOR, Key for IRS measurement , Col , (3- 13) . Synthetic photometry, IR,S data points are averaged in 11 adj acent bins , Wc 
used 2/^m wide bins at A < 20/j,m and 3/^m wide bins at A > 20. The quoted uncertainties include both the instrumental noise and the variation of the SED within a bin. 
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TABLE 6 
Disk properties 



Source ID Tduat [K] Rdust [AU] fdu>t [lO""*] Reduced 



STiLJ oO ( U 


Ooi 1 


/I 9-l-Q 


1^ A -l-n A 

0.4l:iU.41: 


n Q 
u.y 


X 1 LJ J.OUOU 


<^ D4: 


^01 


~0 16 






61it 1 


36it 1 


5 libO 2 


3 4 


rlU 10/ 41:0 




1 S-l- 1 




9 T 


±u 1 41:0 


62it 1 


47i2 


3 8ib0 2 


7 6 








9 1 -i-n 1 


9 s 




1 1 R-t- 1^ 


1 ri-l- 1 


1 ns-t-n nfi 


u.o 


nxj 1 u 


51di3 


75i9 




2 


nx^ ou^fi / 


67di 1 


33it 2 


9 2it0 6 


1 8 


T TT~\ Q01 Qc; 


89it9 


12it3 


n fif^-i-n 1 ^ 

U.DO^U. ±0 


1 4 


X 1 LJ ooUo± 


00^41: 


x1 "^-1-7 


n /I fi-l-fi n7 


n 7 


T TT~^ Qc; 1 1 /I 


Q7-I-Q 

y / iy 


12i2 




n A 




69ib 1 


23it5 


15 2itl 


1 4 


xrxj_/ o\ji?\jo 


58di 1 


45 i 9 


13 4itl 


1 5 


rlU DUO i 1 


41:0^.!; 




1. OiU. 1 


n 9 


HD 113337 


53±1 


55±3 


0.98±0.07 


0.5 


HD 120160 


57±4 


84±14 


0.81±0.20 


0.2 


HD 125451 


63±5 


37±6 


0.18±0.02 


0.2 


HD 127821 


45±1 


66±3 


2.1±0.1 


1.3 


HD 151044 


57±2 


32±2 


0.77±0.05 


0.7 


HD 170773 


43±1 


78±3 


4.8±0.2 


1.3 


HD 192758* 


61±1 


45±9 


5.7±0.3 


3.7 


HD 205674 


54±1 


46±3 


3.7±0.3 


2.4 


HD 206893 


49±1 


49±2 


2.5±0.1 


1.7 


HD 213429 


< 62 


>27 


~0.08 




HD 213617 


55±1 


59±3 


0.96±0.05 


0.7 


HD 221853 


84±1 


23±1 


7.9±0.4 


1.6 



Note. — Disk parameters, listed in this table, come from a model assuming a single narrow dust ring. Disks marked by asterisks can be better 
fitted with a two-component model (see Table [TJ. Col.(l): Identification. Disks discovered in this programme are in boldface. Col. (2): Disk 
temperature. Col. (3): Disk radius. Col. (4): Fractional dust luminosity /dust — Col. (5): Best reduced . 



TABLE 7 
Disk properties 



Warm dust Cold dust 

Source ID Tdu.t [K] Rdu=t [AU] fdust [10^''] Tdu^t [K] Rd„st [AU] fdu.t [lO^*] Reduced 



HD 15115 


179±46 


4±2 


0.38±0.08 


57±1 


42±2 


4.8±0.2 


1.6 


HD 15745 


147±22 


6±2 


3.5±1.6 


81±3 


21±2 


18.9±2.0 


1.1 


HD 16743 


147±24 


8±3 


0.51±0.07 


53±1 


63±4 


3.6±0.3 


1.7 


HD 30447 


159±36 


6±3 


0.66±0.32 


62±2 


38±3 


8.8±0.7 


0.8 


HD 192758 


154±31 


7±3 


0.39±0.09 


56±1 


53±11 


5.4±0.3 


1.8 



Note. — Col.(l): Identification. Col. (2): Temperature of warm dust in the inner ring. Col. (3): Radius of the inner dust ring. Col. (4): 
Fractional dust luminosity of the inner dust ring. Col. (5): Temperature of cold dust in the outer ring. Col. (6): Radius of the outer dust ring. 
Col. (7): Fractional dust luminosity of the outer dust ring. Col. (8): Best reduced x^- 
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TABLE 8 
Age estimates 



Source ID 




log 


Ref. 


Membership 


Rcf. 


Ago [Myr] 


Dating method 


R,ef 


rlU oD f U 


— 4.39 






Coluniba assoc. 


1 


30 


1 




rlJJ ioUbU 












230U±iOU 


6 


1 


HD 15115 


—4.93 






• ' ' 

p Pic mg. 


3,4 


12 


1 




rlJJ iO(4o 








/3 Pic mg. 


1 


12 


1 




rlJJ iDi4o 












10—50 


2,3,5 




HD 17390 


— 5.15 












6 


1 


HD 24636 


-5.41 






Tucana-Horologium assoc. 


1 


30 


1 




HD 25570 


-5.24 






Hyadcs? 


sec Seet. 14.51 


625±50 


1 




HD 30447 








Columba assoc. 


3,4 


30 


1 




rlU oziyo 


— 3.92 






Tucana-Horologium assoc. 


4,5 


30 


1 




HD 33081 












3100+5^00 


6 


1 


HD 35114 


-3.86* 






Columba assoc. 


4 


30 


1 




HD 35841 








Columba assoc. 


3,4 


30 


1 




rlD obyDo 








Octans assoc. 


1 


20 


1 




HD 50571 


— 5.27 


—4.55 




B3 group 


1 


300±120 


1,6 


3 


HD 113337 


— 5.07* 










40±20 


2 




HD 120160 












1300±100 


6 


1 


HD 125451 


-5.10 


-4.37 


3 


Ursa Major mg. 


2 


500±100 


1 




HD 127821 


-5.11 










220±50 


6 


2 


HD 151044 




-5.0 


4 








2,3 




HD 170773 


-4.98* 


-4.39 


1 






200 


6 


3 


HD 192758 








Argus assoc. 


3 


40 


1 




HD 205674 


-5.13 






ABDor mg.? 




70-300 


1,6 


3 


HD 206893 


-4.87* 


-4.47 


1 






2QQ+1000 


6 


1,3 


HD 213429 




-4.83 


2 






2200+1300 


3,6 


1 


HD 213617 












1200±300 


6 


1 


HD 221853 








Local Association 


3 


20-150 


1 





Note. — Col.(l): Identification. Col. (2): Fractional X-ray luminosity based on ROSAT data. Asterisks indicate those objects where the 
correlation between the X-ray source and the star is confirmed by observations with the XMM satellite (XMM-Ncwton slew su rvey Source 
Catalogue, Saxton ct al. 20081) as well. Col. (3): Fractional Call H&K luminosity (logi?^j^). Col. (4): References for log data: 1) IGray et al.| 

1^2006), 2) Gray ct al., 120(51) ; 3) IKing Schuleil Ii2005.) : 4) IWright et"al| 120041) . Col. (5). Membe rship status of the sta r. Col. (6). Refere nces 
for the ident i ficatio n of star as a member o f a spec ific kinematic group in Col. (5): 1) this work; 2) IKin£ et all f200^ ; 3) I Moor et all f2Q0^ : 4) 
iTorres et al.| ( I2008D : 5) iZuckerman &z Songl l |2004bl) . Col. (7). Estimated age of the star and its formal uncertainty. Col. (8). Used age dating 
methods: 1) stellar kinematic group or cluster membership; 2) isochrone fitting; 3) diagnostic of chromospheric activity indicator; 4) diagnostic 
of coronal activity indicator; 5) lithium abundance; 6) literature data. Col. (9). References for the literature data mentioned in Col. (7-8): 1) 
IHolmberg et al.| ( 120091) : 2) Moor et al., 2010b, in prep.; 3) IRhee et al.l II2007I) . 
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Fig. 1. — Offset between the source positions at 24/im and the 2MASS position (after correcting for proper motion due to the time 
difference between the observations). Plus signs indicate stars witli pure photosperic emission, circles show stars exhibiting excess at one 
or more MIPS wavelengths, while targets surrounded by bright extended nebulosity at 24/xm and/or 70/xm are represented by diamonds. 
The distributions of the offsets for stars with and without excess are in agreement within their formal uncertainties. 
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Fig. 2. — Flux ratio of the measured to the predicted flux densitites as a function of the predicted photospheric fluxes for our sample stars 
measured at 24^m with MIPS. For symbols see the caption of Fig. [T] Small panel: the histogram of the flux ratio. The peak at around 
unity can be fitted by a Gaussian with a mean of 1.00 and dispersion of 0.038 (dashed line). We note that the disks around HD 15745, 
HD 35841 and HD 221853 with flux ratios of 9.79, 3.39, and 3.87, respectively, are out of the displayed range. 
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Fig. 3. — Positional offsets between the centroids of point sources detected at and the 2MASS position (after correcting for proper 

motion due to the time difference between the observations) as a function of the signal-to-noise ratio measured at TOfim. For symbols see 
the caption of Fig. [l] Apart from HD 14691, HD 48391 and HD 86146 all of the detected sources exhibit excess emission at 70^m. 
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Fig. 4. — Histogram of the significances of the differences between the measured and predicted photospheric flux densitites, defined as 
{F^o — ^Vo)/o'7o*i for targets measured at 70/im. Small panel; a zoom for the peak at zero. A Gaussian fit to the peaJs provides a mean of 
—0.09 and a dispersion of 0.87, in good agreement with the expectations. 
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Fig. 5. — The ratio of the synthetic IRS 24/im photometry to the MIPS photometry at 24/im as a function of the MIPS flux densities. 
For symbols see the caption of Fig.[T] The mean of the IRS2A/ ratios is 1.02 with a dispersion of 0.06. The dashed line corresponds to 
the IRS2A/F2A ratio of 1. 



Evolution of debris disks around F-type stars 



31 




40 20 -20 -40 8 9 10 15 20 30 

RA offset ["] Wavelength [[im] 

Fig. 6. — Left: 2MASS Ks image of HD 38905 and its surroundings. The position of the two Spitzer/IRS slits arc ovcrplotted. 
A nearby source at a distance of 13.5" is included in the chO slit. At these wavelengths, the two sources are well resolved, 
and a separate spectrum for each source can be extracted. The ch2 slit is nearly perpendicular to the chO slit, but due to 
the longer wavelength, the PSF is wider, thus the nearby source might have a contribution to the ch2 spectrum extracted for 
HD 38905. Right: the solid black line indicates the spectrum extracted at the position of HD 38905, and the black dot is a MIPS 
24 /im photometric point for HD 38905. The black dotted line indicates the spectrum of the nearby source: for A < 14.5 /im, 
it is observed and resolved by IRS, and the black square is also a resolved MIPS 24 /xm photometric point. The straight line 
above 14.5 /ttm is a linear extrapolation. Using the IRS beam profiles and the spectrum of the nearby source, we estimated the 
contribution this source has in the ch2 slit (gray dash-dotted line). The gray dashed line represents the stellar photosphere of 
HD 38905. The black dash dotted line is the sum of the photosphere and the contribution of the nearby source. Our conclusion 
is that the excess emission with respect to the stellar photosphere observed at the position of HD 38905 can be well explained by 
the contamination from the nearby source. The shape of the nearby source's spectrum {Fv ~ A) indicates that it is probably a 
background galaxy (see e.g. Wu ct al. 2009 or Buchanan ct al. 2006). A similar analysis was done for HD 34739, HD 145371 and 
HD 184169, although in those cases, no resolved spectroscopy is available for the nearby sources. Supposing that these sources 
are also background galaxies, a spectral shape of F,^ ^ A is assumed and absolute brightness level was scaled to resolved 24 fim 
MIPS photometry. Our analysis indicates that apart from these four stars, no other targets have suffered contamination by 
nearby sources. 



32 



Moor et al. 



3.0 



00 



2.5 



2.0 



[] 




HD 170773 



□ HD 50571 I 



20 



40 



60 

SNR. 



80 



100 



120 



Fig. 7. — The ratio of the flux density measured in apertures with radii of 18" and 8" as a function of the SNR obtained in the smaller 
aperture. Squares indicate our stars that exhibit excess at 70^m, triangles correspond to debris disks (HD 10647, HD 388 58, HD 48682, 
HD 10 5211, HD 115617, HD 109085, HD 139664, HD 207129) that found to be marginally resolved at this wavelength bv IBrvden eTaTI 
ll200g) . 
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Derived /3 values for some selected disks (see Sect. 14.41 



34 



Moor et al. 



„ 10' 

>. 10"' 

I 10° 

i 10-' 

„ 10' 

>> 10 ' 

« 10 = 

I 10-' 

„ 10' 

>i 

;^ 

>. 10"' 

I 10"° 

I lo-- 

„ 10' 

>> 10' 

^ 10"' 
X 

E 10-' 

„ 10' 

>i 

[3 

5 10"' 

I 10"° 

I 10"= 

„ 10' 

>> 10"' 

^ 10"» 

I 10"= 

„ 10' 

5 10"' 

I 10"° 

E 10-' 

„ 10' 

>■ 10"' 

c ., 

^ 10"' 

I 10-' 

„ 10' 
[3 

5 10"' 

I 10° 

I 10"= 



HD 3670 


HD 15060 


HD 15115 

^^^^^ 


' 'HD 15745 ' - 


' kD 16743 ' - 


' kO 17390 ' " 
""■■■I 

^^^^ 


'hD 24636 ' - 


- ' kD 25570 ' - 


■ ' kD 30447 ' ■ 


■' ' 'hD 32195 ' - 


- ' kD 33081 ' - 


■ ' kD35ll4 ' ■ 


'hD 35841 ' - 


kD 36968 ' - 
...©■•■0 


kD 50571 ' " 


■ ' 'hB 113337 ' - 


- ' kD 120160 ' - 

;---.<^.,..., ; 


■ ' kD 125451 ' ■ 


to 127821 ' - 
^"'^--^ 


kO 161044 ' - 


kD 170773 ' " 


' 'hD 192758 ' - 


- ' kD 205674 ' - 


■ ' kD 206893 ' ■ 


'HD 213429 ' - 


kD 213617 ' - 


kD 221853 ' 



10 100 1000 10 100 1000 10 100 1000 

Wavelength [^m] Wavelength [^m] Wavelength [jim] 



Fig. 9. — Spectral energy distributions for stars exhibiting IR excess in our sample. The different symbols represent the 
folio-wing photometric data - red circles: MIPS, blue squares; ISOPHOT, yello-w triangles: IRAS, blue upside do-wn triangles: 
submillimeter /millimeter observations. The IRS and MIPS SED spectra are displayed -with green lines. The photospheric models 
and the disk models are shown by solid grey lines and dotted black lines, respectively. For HD 15115, HD 15745, HD 16743, 
HD 30447 and HD 192758 the two-component disk models are displayed. 
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Fig. 10. — a) The distribution of Li equivalent widths as a function of V — Kg color indices, b) The distribution of fractional x-ray 
luminosities as a function of V — Kg color indices, c) Color-magnitude diagram. Grey circles represent the known members of the /3Pic 
moving group. Grey squares represent binary stars in the GMD. The black diamond represents BD-|-45° 598, the black triangle shows 
HD 15745. 
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Fig. 12. — H-R diagram for HD 113337B. The effective temperature of the star and its uncertainty are estimated based on spectral 
class information (1 subclass uncertainty is assumed), the absolute magnitude in K-band is computed by assuming that it is located at 
the same distance as the primary component. Mass tracks and isochrones computed by Siess et al. (2000) for the evolution of stars with 
solar metallicity are overplotted. Based on these evolutionary tracks we find the age of the HD 113337B to be 40±20Myr and the mass is 
~0. 13-0.25 Mq. 
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Fig. 13. — /S's-[24] color versus the effective te mperature. Tlie squar e sliows the position of IfD 113337B, the dashed line represents the 
locus of stellar photospheric colors determined bv lGautier et all l|2007l ). 
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Fig. 14. — Histogram of the derived dust temperatures. In those cases where the presence of two separated dust rings are assumed only 
the temperature of the colder ring has been tal<en into account. 
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Fig. 15. — Fractional luminosity of the infrared excess as a function of age. Typical uncertainties in fractional luminosity range from 
0.02 dex to 0.12 dex. Crosses mark disks where the excess emission was measured only in one IR band, thus they have less reliable 
fractional luminosities. Models of de bris disk evolution predict that the decay of fractional luminosity is proportional to t~°' where a 
ranges between 0.3 and 1.0 (see Sect. l5.4ll . For comparison with our data, we plotted the two extremes of the evolutionary models with 
arbitrary normalization (dashed lines). The distribution of the data points, in particular their upper envelope, seems to suggest a decay 
rate halfway between the two extremes. 
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Fig. 16. — Derived dust radii as a function of age. Typical formal uncertainties in the derived disk radii range from 5% to 20%. The 
inset shows a comparison between the cumulative distribution of disks' radii around stars with age <100Myr and stars with age >100Myr 
(disks with lower limits for radius are not included). The inset does not cover any symbols. 
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Fig. 17. — The relative dust masses (/dust ■ ^dust Sect I5.5ll for the disks with ages of 30Myr as a function of the derived radii. 
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Fig. 18. — Histogram of the minimum ^'m,min values that are required for the observed disk to become self stirred in t <C isystc 
disk-to-star mass ratios corresponding to the computed Xm.min are indicated on the top of the graph (Sect. [575)l . 
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